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Basic mechanisms controlling orbital order and orbital fluctuations in transition metal 
oxides are discussed. The lattice driven classical orbital picture, e.g. like in manganites 
LaMnC>3, is contrasted to the quantum behavior of orbitals in frustrated superexchange 
models as realised in pseudocubic titanites ATiOs and vanadates AVO3. In YVO3, the 
lattice and superexchange effects strongly compete - this explains the extreme sensitivity 
of magnetic states to temperature and doping. Lifting the ti g orbital degeneracy by a 
relativistic spin-orbital coupling is considered on example of the layered cobaltates. We find 
that the spin-orbital mixing of low-energy states leads to unusual magnetic correlations in 
a triangular lattice of the C0O2 parent compound. Finally, the magnetism of sodium-rich 
compounds Nai_ a! Co02 is discussed in terms of a spin/orbital polaronic liquid. 

§1. Introduction 

It is now well recognized that the orbital degrees of freedom are an important 
control parameter for physical properties of transition metal oxidesP^® The or- 
bital quantum number specifies the electron density distribution in crystal, hence 
providing a link between magnetism and structure of the chemical bonds3'® When 
the symmetry of the crystal field experienced by a magnetoactive <i-electron is high, 
we encounter the orbital degeneracy problem. In case of a single impurity this 
degeneracy, being an exact symmetry property, cannot be lifted by orbital-lattice 
coupling; instead, a purely electronic degeneracy is replaced by the same degeneracy 
of so-called vibronic states, in which the orbital dynamics and lattice vibrations are 
entangled and cannot be longer separated. However, the orbital degeneracy must 
be lifted in a dense system of magnetic ions. The basic physical mechanisms that 
quench the orbital degrees of freedom in a solid — via the orbital-lattice Jahn- Teller 
(JT) coupling, via the superexchange (SE) interactions between the orbitals, and 
via the relativistic spin-orbit coupling — are well known and have extensively been 
discussed earlier®® see, in particular, a review article.® The purpose of this paper 
is to discuss some recent developments in the field. 

Usually, it is not easy to recognize which mechanism, if any, is dominant in a 
particular material of interest. As a rule, one expects strong JT interactions for 
the orbitals of e g -symmetry as they are directed towards the ligands. On the other 
hand, Ti, V, Ru, etc. ions with t^g orbital degeneracy are regarded as "weak Jahn- 
Teller" ions, and the other mechanisms listed above might be equally or even more 
important in compounds based on these ions. An empirical indication for the strong 
JT case is a spin-orbital separation in a sense of very different (magnetic T m and 
structural T s t r ) transition temperatures and large lattice distortions observed. Or- 
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bitals are (self)trapped by these coherent lattice distortions, and (i-shell quadrupole 
moments are ordered regardless to whether spins are ordered or not. While the 
other two mechanisms are different in this respect and better characterized, loosely 
speaking, by spin-orbital confinement indicated by T m ~ T str . Indeed, strong cou- 
pling between spin and orbital orderings/fluctuations is an intrinsic feature of the 
superexchange and spin-orbit interactions "by construction". The superexchange 
mechanism becomes increasingly effective near the Mott metal-insulator transition, 
because the intensity of virtual charge fluctuations (which are ultimately responsible 
for the exchange interactions) is large in small charge-gap materials. An increased 
virtual kinetic energy of electrons near the Mott transition — in other words, the 
proximity to the metallic state — can make the electronic exchange more efficient 
in lifting the orbital degeneracy than the electron-lattice coupling.^ ^ 

As the nature of interactions behind the above mechanisms are different, they 
usually prefer different orbital states and compete. A quantitative description of 
the orbital states is thus difficult in general but possible in some limiting cases. In 
particular, if the ground state orbital polarization is mainly due to strong lattice 
distortions, one can just "guess" it from the symmetry point of view, as in standart 
crystal field theory. Low-energy orbital fluctuations are suppressed in this case, 
and it is sufficient to work with spin-only Hamiltonians operating within the lowest 
classical orbital state. As an example of such a simple case, we consider in Section 
2 the magnetic and optical properties of LaMn03. 

The limit of strong (superexchange and relativistic spin-orbital) coupling is more 
involved theoretically, as one has to start in this case with a Hamiltonian which 
operates in the full spin-orbital Hilbert space, and derive the ground state spin- 
orbital wavefunction by optimizing the intersite correlations. It turns out that this 
job cannot be completed on a simple classical level; one realizes soon that the spin- 
orbital Hamiltonians possess a large number of classical configurations with the same 
energy. Therefore, theory must include quantum fluctuations in order to lift the 
orbital degeneracy and to select the ground state, in which the orbitals and spins 
can hardly be separated — a different situaion to compare with the above case. The 
origin of such complications lies in the spatial anisotropy of orbital wavefunctions. 
This leads to a specific, non-Heisenberg form of the orbital interactions which are 
frustrated on high-symmetry lattices containing several equivalent bond directions. 
Physically, orbital exchange interactions on different bonds require the population of 
different orbital states and hence compete. This results in an infinite degeneracy of 
the classical ground states. A substantial part of this paper is devoted to illustrate 
how the orbital frustration is resolved by quantum effects in spin-orbital models, and 
to discuss results in the context of titanites and vanadates with perovskite structure 
{Section 3). We will also demonstrate a crucial difference in the low-energy behavior 
of superexchange models for e g and t2 g orbitals. 

In some cases, the competition between the superexchange and orbital-lattice 
interactions result in a rich phase diagram, including mutually exclusive states sep- 
arated by first order transitions. A nice example of this is YVO3 discussed in Sec- 
tion 4- In particular, we show how a competition between the three most relevant 
spin/orbital electronic phases is sensitively controlled in YVO3 by temperature or 
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small amounts of hole-doping. 

The last part of the paper, Section 5, discusses the role of the relativistic spin- 
orbit coupling X so , which effectively reduces the orbital degeneracy already on the 
single-ion level. This coupling mixes-up the spin and orbital quantum numbers in the 
ground state wavefunction, such that the magnetic ordering implies both the orbital 
and spin ordering at once. The spin-orbit coupling mechanism might be essential in 
insulating compounds of late-3<i ions with large X so and also in 4d-ruthenates. We 
focus here on two-dimensional cobaltates with a triangular lattice, and present a the- 
ory which shows that very unusual magnetic states can be stabilized in C0O2 planes 
by spin-orbit coupling. We also discuss in that section, how the orbital degeneracy 
and well known spin-state flexibility of cobalt ions lead to a polaron-liquid picture 
for the sodium-rich compounds Nai-zCoC^, and explain their magnetic properties. 

Our intention is to put more emphasis on a comparative discussion of differ- 
ent mechanisms of lifting the orbital degeneracy, by considering their predictions 
in the context of particular compounds. Apart from discussing previous work, the 
manuscript presents many original results which either develop known ideas or are 
completely new. These new results and predictions made may help to resolve con- 
troversial issues and to discriminate between the different models and viewpoints. 

§2. Lifting the orbital degeneracy by lattice distortions 

Let us consider the Mott insulators with composition ABO3, where A and B 
sites accommodate the rare-earth and transition metal ions, respectively. They crys- 
tallize in a distorted perovskite structure P®' 1 Despite a very simple, nearly cubic 
lattice formed by magnetic ions, a variety of spin structures are observed in these 
compounds: An isotropic G-type antiferromagnetism (AF) in LaTiC-3, isotropic fer- 
romagnetism (F) in YTi03, and also anisotropic magnetic states such as C-type AF 
in LaV03 and ^4-type AF in LaMn03. 101 The richness of the spin-orbital states, 
realized in these compounds with a similar lattice structure, already indicates that 
different mechanisms lifting the orbital degeneracy might be at work, depending on 
the type of orbitals, the spin value, and on the closeness to the Mott transition. 

Considering the lattice distortions in perovskites, one should distinguish distor- 
tions of two different origins: (A) The first one is due to ionic-size mismatch effects, 
that generate cooperative rotations and also some distortions of octahedra in order 
to fulfill the close-packing conditions within a perovskite structure. These are the 
"extrinsic" deviations from cubic symmetry, in a sense that they are not triggered 
by orbitals themselves and are present even in perovskites having no orbital degen- 
eracy at all, e.g. in LaA103 or LaFeC>3. The orbitals may split and polarize under 
the extrinsic deformations, but they play essentially the role of spectators, and to 
speak of "orbital ordering" in a sense of cooperative phenomenon (such as "spin 
ordering") would therefore be misleading in this case. (B) Secondly, cooperative 
JT-distortions, which are generated by orbital-lattice coupling itself at the orbital 
ordering temperature (seen by "eyes" as a structural transition). 

Usually, these two contributions are superimposed on each other. Sometimes, it 
is not easy to identify which one is dominant. The temperature dependence of the 
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distortion is helpful in this context. Manganites show an orbital order-disorder phase 
transition at high temperatures at about 800 K, below which a large, of the order of 
15%, distortion of octahedra sets in. This indicates the dominant role of a cooperative 
JT physics, which is natural for e g orbital systems with strong coupling between 
the oxygen vibrations and e g quadrupole. In contrast, no cooperative structural 
phase transition has thus far been observed in titanates, suggesting that small lattice 
distortions present in these compounds are mostly due to the ionic-size mismatch 
effects. This seems also plausible, as i2g-orbital JT coupling is weak. 

Whatever the origin, the lattice distortions generate low (noncubic) symmetry 
components in the crystal field potential, which split the initially degenerate orbital 
levels. While structure of the lowest, occupied crystal-field level can be determined 
without much problem — simply by symmetry considerations, the level splittings 
are very sensitive to the value of (i) distortions and (ii) the orbital-lattice coupling 
(both factors being much smaller in tig orbital compounds). If the splittings are large 
enough to localize electrons in the lowest crystal-field level and to suppress the inter- 
site orbital fluctuations, classical treatment of orbitals is justified. Accordingly, the 
e 3 -electron density is determined by a function parameterized via the site-dependent 
classical variables af. 

ifti = cosaj|3z 2 — r 2 > + sinaj|a; 2 — y 2 > , (2T) 

while the occupied ti g orbital can be expressed via the two angles: 

ipi = cos cti\xy > + sin«j cos (3i\yz > + sinoj sin/3j|zx > . (2 - 2) 

A knowledge of these functions allows one to express a various experimental observ- 
ables such as spin interactions, optical absorption and Raman scattering intensities, 
etc. in terms of a few "orbital angles". [In general, classical orbital states may 
include also a complex wave-functions, but they are excluded in the case of strong 
lattice distortions^^]. These angles can thus be determined from the experimental 
data, and then compared with those suggested by a crystal-field theory. Such a clas- 
sical approach is a great simplification of the " orbital physics" , and it has widely and 
successfully been used in the past. Concerning the orbital excitations in this picture, 
they can basically be regarded as a nearly localized transitions between the crystal- 
field levels. We demonstrate below how nicely this canonical "orbital-angle" scenario 
does work in manganites, and discuss its shortcomings in titanites and vanadates. 

2.1. e g orbitals: The case of LaMnO^ 

Below a cooperative JT phase transition at ~800 K^> T/v, a two-sublattice 
orbital order sets in. Suggested by C-type arrangement of octahedron elongation 
(staggered within ab planes and repeated along c direction), the lowest crystal field 
level can be parameterized via the orbital angle ffi as follows: 

|± >= cos ^ |3-z 2 - r 2 > ±sin^|:E 2 - y 2 > . (2-3) 

In this state, the e g electron distribution is spatially asymmetric, which leads to 
strong anisotropy in spin exchange couplings and optical transition amplitudes. 
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Thus, the orbital angle 9 can be determined from a related experiments and com- 
pared with 6 ~ 108^ suggested by structural data. 

Spin interactions. — For manganites, there are several intersite virtual transitions 
dfdj — ► dfdj, each contributing to the spin-orbital super exchange. The correspond- 
ing energies can be parameterized via the Coulomb repulsion U = A + 4B + 3C for 
electrons residing in the same e g orbital (Un a ^n a i), the Hund's integral Jh = 45+C 
between e g spins in different orbitals (— 2Jhs q - sp), the Hund's integral J' H = 2B + C 
between the e g spin and the t 2g -coYe spin (—2J' H St • s e ), and, finally, the Jahn- Teller 
splitting Ajt between different e 9 -orbitals: Ajxinp — n a )/2. (A, B, C are the Racah 
parameters, see RefPJ for details). For the e g electron hopping dfdj — > dfd^j, we 
obtain the following five different transition energies: 

E 1 = U + A JT - 3J H , (2-4) 
E 2 = U + Aj T - 3J H + 5J' H , 



E 3 = U + Aj T + 3J' H - yfA 2 JT + J2 , 

E 4 = U + Aj T + 3J' H - J H , 

E 5 = U + Aj T + 3J' H + ^Aj T + J 2 H . 

There are also (weaker) transitions associated with t 2g electron hoppings, which 
provide an isotropic AF coupling between the spins J t (Si -Sj), with J t = t' 2 /2E t and 
S = 2 of Mn 3+ . Here, E% = U + Jh + 2J' H is an average excitation energy for the 
intersite t 2g — electron hoppings, and t' ~ t/3 follows from the Slater-Koster relation. 

From a fit to the optical conductivity in LaMnOs, the energies E n have been 
obtainedPJ Then, it follows from Ens. ([2^ that U ~ 3.4 eV, J' H ~ 0.5 eV, J H ~ 
0.7 eV and Ajt ~ 0.7 eV (in the present notations for U, J' H , Jh given above; Refill 
used instead U = A—2B+3C and Jh = 2B+C). The Hund's integrals are somewhat 
reduced from the atomic values J' H = 0.65 eV and Jh = 0.89 eVp) The e 9 -orbital 
splitting is substantial, suggesting that orbitals are indeed strongly polarized by 
static JT distortions, and justifying a crystal-field approach below 800 K. A large 
JT binding energy (= Ajt/ '4) also indicates that dynamical distortions of octahedra 
are well present above 800 K, thus a structural transition is of the order-disorder type 
for these distortions and e g — quadrupole moments. This is in full accord with a view 
expressed in RefP Note, however, that JT energy scales are smaller than those set 
by correlations, thus LaMnOs h as to be regarded as a typical Mott insulator. 

The superexchange Hamiltonian H consists of several terms, Hn \ originating 
from virtual hoppings with energies E n Q2-4D in the intermediate state (7 denotes 
the bond directions a, b, c). Following the derivation of Ref.pl we obtain: 

H v = & - • ^ + 6)(1 - 4r ^ ,)(7) + ( i + si )( ^ • §j - 4)(1 - 4r ^ ,)(7) 
+ ( si + 8^ )( ^ ' ^ " 4)(1 " 2Ti)(7)(1 " 2Tj)(7)] ' (2 ' 5) 



where the e 5 -pseudospins are defined as 2t- c ^ = of, 2t^ = cos (fiaf + sin 4>af, 
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2t\ = cos (fxrf — sin (fiof with tfi = 2tt/3. Here, o z and a x are the Pauli matrices, and 
a z =l (-1) corresponds to a \x 2 — y 2 ) (\3z 2 — r 2 }) state. Note that no u^-component is 
involved in operators; this relfects the fact that the orbital angular momentum 
is fully quenched within the e g doublet. Physically, a pseudospin t^> —structure of 
H2-5|) reflects the dependence of spin interactions on the orbitals that are occupied, 
expressing thereby a famous Goodenough-Kanamori rules in a compressed formal 
wayP 

In a classical orbital state (|2-.3|) . operators a a can be regarded as a numbers 
a z = — cos0 and a x = ± sin#, thus we may set (1 — AriTj)^ = sin 2 9, (1 — 4TiTj) (afe ) = 
(3/4 + sin 2 9) etc. in Eq.(j23J. At this point, the spin-orbital model (f2^|) "collapses" 
into the Heisenberg Hamiltonian, J^(9)(Si ■ Sj), and only "memory" of orbitals that 
is left is in a possible anisotropy of spin-exchange constants J a ,b,c'- 

J " = t K - i + i + i> <i + 8in2 ") + 1 <i + (I - cos 9 ' 2 ] + J - 

Note that the £^i-term has a negative (ferromagnetic) sign. This corresponds to the 
high-spin intersite transition with the lowest energy E%, which shows up in optical 
absorption spectra as a separate band near 2 eVp 2 ^ From the spectral weight of 
this line, one can determine the value of t, see below. All the other terms come 
from low-spin transitions. The orbital angle 9 controls the competition between a 
ferromagnetic E\— term and AF E2 r .-,E^— contributions in a bond-selective way. 

Optical intensities. — Charge transitions d\d^ — > dfd^ are optically active. The 
spectral shape of corresponding lines are controlled by band motion of the excited 
doublons (holes) in the upper (lower) Hubbard bands and by their excitonic binding 
effects. The intensity of each line at E n is determined by a virtual kinetic energy, 
and thus, according to the optical sum rule, can be expressed via the expectation 
values of superexchange terms Hn\ij) for each bond direction 7. 

The optical intensity data is often quantified via the effective carrier number, 

(7 ) 2m wo f°° (y) 1 v , (n „s 

N eff,n = -^rj Q W(")du>, (2-7) 

where mo is the free electron mass, and v$ = a|] is the volume per magnetic ion. 
Via the optical sum rule applied to a given transition with n = 1, 5, the value of 
^e//n can ' 3e ex P resse d as follows?^ 

^]l = ^K^ = -^(2H^m- (2-8) 

Here, is the kinetic energy gain associated with a given virtual transition n for 
a bond (ij) along axis 7. The second equality in this expression states that if„ 7 , 
hence the intensity of a given optical transition, is controlled by expectation value of 
the corresponding term in spin-orbital superexchange interaction ()2-5|) , Via the 
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operators , each optical transition E n obtains its own dependence on the orbital 
angle 9. Thus, Ea. (|2-8j) forms a basis for the study of orbital states by means of 
optical spectroscopy, in addition to the magnetic data. Using a full set of the optical 
and magnetic data, it becomes possible to quantify more reliably the values of t, 
U, Jh in crystal, and to determine the SE-energy scale, as it has been proposed in 
Ref.f^ and worked out specifically for LaVO^ and LaMnOsPJ 

Physically, spin and orbital correlations determine the optical intensities of dif- 
ferent transitions E n via the selection rules, which are implicit in Ea. (|2-5|) . For 
instance, the intensity of the high-spin transitions obtained from the Ei-term in 
Eq.(H2J) reads: 

K i ] = T^r <Si-Sj + 6 >® sin 2 9, (2-9) 
K i b) = TtSt < ^ ■ S 3 + 6 > W (3/4 + sin 2 9). 

At T < T N ~ 140K, < Si-Sj >( afe )^ 4 and < Si-Sj -4 for the A-type classical 

spin state, while < Si ■ Sj >=^ OatT> T/v- Thus, both spin and orbital correlations 
induce a strong polarization and temperature dependence of the optical conductivity. 
Note, that the e s -hopping integral t determines the overall intensity; we find that 
t ~ 0.4 eV fits well the observed optical weig htsPJ This gives 4t 2 /U ~ 190 meV 
(larger than in cuprates, since t stands here for the 3z 2 — r 2 orbitals and thus is 
larger than a planar x 2 — y 2 orbital-transfer). 

A quantitative predictions of the above " orbital-angle" theory for spin-exchange 
constants and optical intensities, expressed in Eas. (|2-6j) and (|2-9j) solely in terms of 
a classical variable 9, have nicely been confirmed in RefP** It is found that the angle 
9 = 102° well describes the optical anisotropy, and gives the exchange integrals J c = 
1.2 meV, J a b = —1.6 meV in perfect agreement with magnon dataPI Given that 
"the lattice angle" 108° has been estimated from the octahedron distortions alone, 
and thus may slightly change in reality by GdFeO"3-type distortions and exchange 
interactions, one may speak of a quantitative description of the entire set of data 
in a self-contained manner (everything is taken from the experiment). This implies 
the dominant role of (JT) lattice distortions in lifting the orbital degeneracy in 
manganites as expected. Of course, the situation changes if one injects a mobile 
holes, or drives a system closer to the Mott transition. The orbital order is indeed 
suppressed in LaMnOs under high pressure,^ in spite the insulating state is still 
retained. Pressure reduces the ionic-size mismatch effects, and, more importantly, 
decreases the charge gap and thus enhances the kinetic energy. The latter implies an 
increased role of superexchange interactions which, as discussed later, are strongly 
frustrated on cubic lattices; consequently, a classical orbital order is degraded. In 
addition, a weak e g — density modulation like in case of CaFeOjP may also contribute 
to the orbital quenching near the metal-insulator transition. 

2.2. tig orbitals: Titanates 

Depending on the A-site ion radius, two types of magnetic orderings have thus 
far been identified in ATiC>3 compounds: (I) G-type AF state as in LaTi03, which 
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changes to (II) ferromagnetic ordering for the smaller size A-ions as in YTi03, with 
strong suppression of transition temperatures in between.^ Surprisingly enough, 
both AF and F states are isotropic in a sense of equal exchange couplings in all three 
crystallographic directions EDMS Such a robust isotropy of spin interactions (despite 
the fact that the c-direction bonds are structurally not equivalent to those in the ab 
plane) , and a kind of " soft" (instead of strong first order) transition between such a 
incompatible spin states is quite unusual. 

Another unique feature of titanites is that, different from many oxides, no coop- 
erative structural phase transition is observed so far in LaTiC>3 and YTi03 . (Except 
for a sizable magnetostriction in LaTi03, 21 ' indicating the presence of low-energy 
orbital degeneracy). This is very unusual because the JT physics could apparently 
be anticipated for Ti 3+ ions in nearly octahedral environment. One way of thinking 
of this paradox is th&y^Mffl^lMiMl the titanites are located close to the Mott 
transition, and the enhanced virtual charge fluctuations, represented usually in the 
form of spin-orbital superexchange interactions, frustrate and suppress expected JT- 
orbital order (see next section). Yet another explanation is possible,^ which 
is based on local crystal field picture. In this view, the orbital state is controlled 
by lattice, instead of superexchange interactions. Namely, it is thought that extrin- 
sic lattice distortions caused by ionic-size mismatch effects (so-called GdFeC>3-type 
distortions) remove the orbital degeneracy for all temperatures, thus orbital dynam- 
ics is suppressed and no Jahn- Teller instability can develop, either. One caveat in 
this reasoning, though, is that a cooperative orbital transitions are not prevented 
by GdFe03-distortions (of the similar size) in other ti g compounds, e.g. in pseu- 
docubic vanadates. Thus, the titanites seem to be unique and require other than the 
" GdFeOs-explanation" for the orbital quenching. Nonetheless, let us consider now in 
more detail the predictions of such a crystal-field scenario, that is, the "orbital-angle" 
picture expressed in Ea. (|2-2[) . 

The perovskite structure is rather "tolerant" and can accommodate a variation 
of the A-site ionic radius by cooperative rotations of octahedra. This is accompanied 
by a shift of A-cations such that the distances A-0 and A-Ti vary somewhat. Also 
the oxygen octahedra get distorted, but their distortion (~ 3%) is far less than 
that of the oxygen coordination around the A-cation. The nonequivalence of the 
Ti-A distances and weak distortions of octahedra induce then a noncubic crystal 
field for the Ti-ion. In LaTi03, the octahedron is nearly perfect, but the Ti-La 
distance along one of the body diagonals is shorter,^ suggesting the lowest crystal- 
field level of ~ \xz + yz ± xy)/\^3 symmetry, as has been confirmed in Refs.^^ 
by explicit crystal-field calculations. This state describes an orbital elongated along 
the short La-Ti-La diagonal, and the sign (±) accounts for the fact that the short 
diagonal direction alternates along the c-axis (reflecting a mirror plane present in 
crystal structure). Thus, lattice distortions impose the orbital pattern of A- type 
structure. In YTi03, the A-site shifts and concomitant distortions of octahedra, 
induced^ by ionic-size mismatch effects including A-0 and A-Ti covalency, are 
stronger. A crystal-field^*^ and band structure 1 ^ calculations predict that 
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these distortions stabilize a four-sublattice pattern of planar orbitals 

ipi,3 ~ \xz ± xy) /\/2 and ip2,4 ~ \u z ± xy)j\pl . (2-10) 

Note that xy orbitals are more populated (with two times larger filling-factor when 
averaged over sites). This state would therefore lead to a sizable anisotropy of the 
electronic properties. We recall that all the above orbital patterns are nearly temper- 
ature independent, since no cooperative JT structural transition (like in manganites) 
is involved here. 

Earlier neutron diffraction^ and NMRP^ experiments gave support to the 
crystal-field ground states described above; however, this has recently been recon- 
sidered in Refp2' Namely, NMR-spectra in YT1O3 show that the t2 9 -quadrupole 
moment is in fact markedly reduced from the predictions of crystal-field theory, and 
this has been attributed to quantum orbital fluctuations.^ On the other hand, this 
maybe not easy to reconcile with the earlier "no orbital- fluctuations" interpretation 
of NMR spectra in less distorted LaTi03p^ — thus no final conclusion can be made 
at present from NMR data aloneP^ 

For the first-principle study of effects of lattice distortions on the orbital states 
in titanites, we refer to the recent worlP^ in which: the ground-state orbital po- 
larization, (crystal-field) energy level structure, and spin exchange interactions have 
been calculated within the LDA+DMFT framework. In this study, which attempts 
to combine a quantum chemistry with strong correlations, a nearly complete con- 
densation of electrons in a particular single orbital (an elongated, of the a\ g type in 
LaTi03 and a planar one in YTi03, as just discussed) is found, in accord with a sim- 
pler crystal-field calculations of RefsP^*^' Other (empty) orbital levels are found 
to be located at energies as high as ~ 200 meV; this qualitatively agrees again with 
Refs.p^ 1 ^ but deviates from the predictions of yet another first-principle study of 
Ref.ED The ca i cu i a tions of RefP^ give correct spin ordering patterns for both com- 
pounds (for a first time, to our knowledge; — usually, it is difficult to get them both 
right simultaneously, see RefP"). However, the calculated spin exchange constants: 
J c = -0.5 meV, J ab = -4.0 meV for YTi0 3 , and J c = 5.0 meV, J ab = 3.2 meV 
for LaTiOs, — as well as their anisotropy ratio J c / J a b are quite different from 
the observed values: J c ~ J a b ~ —3 meV (YTiOj^J) and J c ~ J a b ~ 15.5 meV 
(LaTi03^). Large anisotropy of the calculated spin couplings J c /J a b ~ 0.1 in 
YTiOs is particularly disturbing; the same trend is found in Ref.: 33 ' J c /Jab ~ 0.3. 
The origin of anisotropy is — put it simple way — that the lattice distortions make 
the c axis structurally different from the two others, and this translates into J c 7^ J a b 
within a crystal-field picture. In principle, one can achieve J c j J a b = 1 by tuning the 
orbital angles and Jh by "hand", but, as demonstrated in Ref.,™ this is highly 
sensitive procedure and seems more like as an "accident" rather than explanation. 
Given that the degree of lattice distortions in LaTi03 and YTiOs are differentP 
- as reflected clearly in their very different crystal-field states — the observation 
Jc/Jab — 1 in both of them is enigmatic, and the results of first-principle calculations 
make the case complete. 

A way to overcome this problem is — as proposed in RefsP^ 1 ^ on physical 
grounds — to relax the orbital polarization: this would make the spin distribution 
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more isotropic. In other words, we abandon the fixed orbital states like (|2-10)) . and 
let the orbital angles to fluctuate, as they do in case of a dynamical JT problem. 
However, a dynamical aspects of JT physics as well as the intersite superexchange 
fluctuations arc suppressed by a static crystal-field splittings, which, according to 
Ref.p^J are large: the ti g level structure reads as [0, 200, 330] meV in YTiC-3, and 
[0, 140, 200] meV in LaTi0 3 . 

Such a strong noncubic splitting of ti g triplet, exceeding 10% of the cubic 
lODg ~ 1.5 — 2.0 eV is quite surprising, and, in our view, is overestimated in the 
calculations. Indeed, e g orbital splitting in manganites has been found ~ 0.7 eV 
(see above) which is about 30% of WDq; this sounds reasonable because of (i) very 
large octahedron distortions (~ 15%) that (ii) strongly couple to the e g -quadrupole. 
Both these factors are much smaller in titanates: (i) octahedron distortions are 
about 3% only, while effects of further neighbors should be screened quite well due 
to a proximity of metallicity; (ii) coupling to the ^-quadrupole is much less than 
that for e g — otherwise a cooperative JT transition would take place followed by a 
strong distortions like in manganites. Putting these together, we "estimate" that a 
t2 g splitting should be at least by an order of value less than that seen for e g level 
in manganites, — that is, well below 100 meV. This would make it possible that 
orbitals start talking to spins and fluctuate as suggested in RefP 

In general, a quantitative calculation of the level splittings in a solid is a rather 
sensitive matter (nontrivial even for well localized /-electrons) . Concerning a numer- 
ical calculations, we tend to speculate that a notion of "orbital splitting" (as used in 
the models) might be not that well defined in first-principle calculations, which are 
designed to obtain the ground state properties (but not excitation spectra). This 
might be a particularly delicate matter in a strongly correlated systems, where effec- 
tive orbital splittings may have strong t/U and frequency dependences.^ As far as 
a simple crystal-field calculations are concerned, they are useful for symmetry analy- 
ses but, quantitatively, often strongly deviate from what is observed. As an example 
of this sort, we refer to the case of Ti 3+ ions in a AI2O3 host, where two optical 
transitions within a t2 g multiplet have directly been observed.^ The level splitting 
~ 13 meV found is astonishingly small, given a sizable (more than 5%) distortion 
of octahedron. A strong reduction of the level splitting from a static crystal-field 
model has been attributed to the orbital dynamics (due to the dynamical JT effect, 
in that case). 

As a possible test of a crystal-field calculations in perovskites, it should be very 
instructive to perform a similar work for LaVC>3 and YVO3 at room temperature. 
The point is that there is a cooperative, of second order orbital phase transition in 
vanadates at low-temperature (below 200 K). To be consistent with this observation, 
the level-spacings induced by GdFe03-type etc. distortions (of the similar size as in 
titanites) must come out very small. 

Whatever the scenario, the orbital excitations obtain a characteristic energy scale 
set by the most relevant interactions. However, the predictions of a local crystal- 
field theory and a cooperative SE interactions for the momentum and polarization 
dependences of orbital response are radically different. Such a test case has recently 
been provided by the Raman light scattering from orbital excitations in LaTiC>3 and 
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YTiO-?, see RefP In short, it was found that the polarization rules, dictated by 
a crystal-field theory, are in complete disagreement with those observed. Instead, 
the intensity of the Raman scattering obeys the cubic symmetry in both compounds 
when the sample is rotated. Altogether, magnon and Raman data reveal a universal 
(cubic) symmetry of the spin and orbital responses in titanites, which is seemingly 
not sensitive to the differences in their lattice distortions. Such a universality is hard 
to explain in terms of a crystal-field picture based — by its very meaning — on the 
deviations away cubic symmetry. 

Moreover, a picture based on a Heisenberg spins residing on a fully polarized 
single orbital has a problem in explaining the anomalous spin reduction in LaTi03 P 
The magnetic moment measured is as small as M ~ 0.5/Us, so the deviation from 
the nearest-neighbor 3D Heisenberg model value Mjj — is unusually large: 

8M/M H = (M H - M)/M H ~ 40%. As a first attempt to cure this problem, one 
may notice that the Heisenberg spin system is itself derived from the Hubbard model 
by integrating out virtual charge fluctuations (empty and doubly occupied states). 
Therefore, the amplitude of the physical magnetic moment is reduced from that 
of the low-energy Heisenberg model via M p h = Ms(l — — rid)- However, this 
reduction is in fact very small, as one can see from the following argument. Let us 
discard for a moment the orbital fluctuations, and consider a single orbital model. By 
second-order perturbation theory, densities of the virtual holons and doublons rid, 
generated by electron hoppings in a Mott insulator, can be estimated as nh = Ud — 
z(jj) 2 , where z is the nearest-neighbor coordination number. Thus, the associated 
moment reduction is SM/Mh ~ ^(tt) 2 - Even near the metal-insulator transition, 
that is for U ~ 2zt, this gives in 3D cubic lattice 8M/M H ~ JL ~ 8% only. (We 
note that this simple consideration is supported by 2D Hubbard model calculations 
of Ref.P the deviation of the staggered moment from 2D Heisenberg limit was 
found SM/Mh — 12% at U ~ 8t, in perfect agreement with the above l/2z— scaling.) 
This implies that the anomalous moment reduction in LaTi03 requires some extra 
physics not contained in a single-orbital Hubbard or Heisenberg models. (In fact, 
an "extra" moment reduction is always present whenever orbital fluctuations are 
suspected: LaTiC>3, LaV03 and C-phase of YVO3). Enhanced spin fluctuations, a 
spatial isotropy of the spin exchange integrals, and similar cubic symmetry found in 
the Raman scattering from orbital fluctuations strongly suggest multiband effects, 
that is, the presence of orbital quantum dynamics in the ground state of titanites. 

2.3. t2g orbitals: Vanadates 

In vanadium oxides AVO3 one is dealing with a t\ g configuration, obtained by a 
removal of an electron from the t% g shell which is orbitally-nondegenerate (because 
of the Hund's rule). From a viewpoint of the JT physics, this represents a hole- 
analogy of titanites: One hole in the t2 g manifold instead of one electron in ATi03. 
Similar to the titanites, they crystallize in a perovskite structure, with GdFe03-type 
distortions increasing from La- towards Y-based compounds, as usual. At lower tem- 
peratures, vanadates undergo a cooperative structural transition (of second-order, 
typically). This very fact indicates the presence of unquenched low-energy orbital 
degeneracy, suggesting that underlying GdFeC>3-type distortions are indeed not suf- 
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ficient to remove it. The structural transition temperature T s t r is nearly confined 
to the magnetic one T/v ~ 140 K in LaVO"3 — the relation T s t r ~ T/v holds even 
in a hole-doped samples,^ — while T s t r ~ 200 K in YVO3 is quite separated from 
T/v ~ 120 K. It is accepted that the xy orbital is static below T str and accommodates 
one of the two magnetic electrons. An empirical correlation between T s t r and the 
A-site ionic size suggests that GdFeC>3-type distortions "help" this stabilization. 

However, a mechanism lifting the degeneracy of xz/yz doublet is controversial: 
Is it also of lattice origin, or controlled by the superexchange? Based on distortions 
of octahedra (albeit as small as in YTiO"3), many researchers draw a pattern of stag- 
gered orbitals polarized by JT interactions. This way of thinking is conceptually 
the same as in manganites (the only difference is the energy scales): a coopera- 
tive, three-dimensional JT-ordering of xz/yz doublets. However, this leaves open 
the questions: Why is the JT mechanism so effective for the ^-hole of vanadates, 
while being apparently innocent (no structural transition) in titanites with one t<i g - 
electron? Why is the T str , if driven by JT physics, so closely correlated with T/v 
in LaV03? Motivated by these basic questions, we proposed a while ago a different 
viewpH* which attributes the difference between vanadates and titanites to the dif- 
ferent spin values: classical S = 1 versus more quantum S = 1/2. Being not much 
important for JT-physics, the spin-value is of key importance for the superexchange 
mechanism of lifting the orbital deg eneracy. An apparent correlation between T s ^ r 
and Tjy in LaVC-3, the origin of C-type spin pattern (different from titanites), etc., 
all find a coherent explanation within the model of RefP^ This theory predicts that 
the xy orbital becomes indeed classical (as commonly believed) below T str , but, dif- 
ferent from the JT scenario, xz/yz doublets keep fluctuating and their degeneracy 
is lifted mainly due to the formation of quasi-lD orbital chains with Heisenberg- 
like dynamics. Concerning the separation of T s j r from T/v in YVO3, we think it is 
due to an increased tendency for the xy orbital-selection by GdFe03-type distor- 
tions; this helps in fact the formation of xz/yz doublets already above Tv- Below 
T s t r , a short-range spin correlations and orbital dynamics is of quasi ID nature, and 
xz/yz doublet polarization is small. Essentially, the xz/yz sector remains almost 
disordered for both quantum and entropy reasons EHJEDE) T n our v iew, a complete 
classical order in the xz/yz sector sets in only below a second structural transition 
at T* tr ~ 77 K, determined by a competition between ~1D spin-orbital physics and 
GdFe03-type distortions, which prefer a different ground state (more to come in 
Section 4). 

Apart from a neutron scattering experiments^ challenging a classical JT pic- 
ture for vanadates, we would like to quote here a recent paperj^ which observes that 
the vanadates become "transparent" for the thermal phonons only below a second 
transition at T* tr (if present), in a so-called low-temperature G-phase, where we in- 
deed expect that everything should be "normal" (described by a static orbital-angle 
physics). Enhanced phonon scattering on xz/yz fluctuations, which suddenly disap- 
pears below T* tr , is very natural within the superexchange model. While it would be 
hard to understand this, if the orbital states both above and below T* tr are classically 
ordered via the JT mechanism. Obviously, the thermal conductivity measurements 
are highly desired in titanites, in order to see as whether the phonon scattering 
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remains unquenched down to low temperatures as expected from a superexchange 
picture, or will be like in manganites with a static orbitals. 

To summarize our present view on the role of orbital-lattice mechanism in per- 
ovskites: (i) The JT physics dominates in manganites, with a secondary role of SE 
interactions and extrinsic GdFeC>3-type distortions; (tt) In titanites and vanadates, 
the "orbital-angle" description is insufficient. It seems that orbital-lattice coupling 
plays a secondary, but important role providing a phase-selection out of a manifold 
of nearly degenerate many-body quantum states, offered by superexchange interac- 
tions. In vanadates with smaller A-site ions, though, a classical G phase, favored 
by GdFeC>3-type distortions, takes over in the ground state, but the quantum spin- 
orbital states are restored again at finite temperature due to their larger spin-orbital 
entropyp^''^ — a somewhat curious cooperation of quantum and thermal effects. 

§3. Lifting the orbital degeneracy by spin-orbital superexchange 

While the kinetic energy of electrons is represented in metals by the hopping 
t-Hamiltonian, it takes a form of spin-orbital superexchange in the Mott insula- 
tor. The superexchange interactions are obtained by eliminating the virtual dou- 
blon/holon states, a procedure which is justified as far as t/U is small, and the 
resulting SE-scale 4t 2 /U is smaller than t itself. Near the Mott transition, a longer- 
range coupling and retardation effects, caused by a softening and unbinding of dou- 
blon/holon excitations are expected, and separation of the spin-orbital modes from 
an emergent low-energy charge and fermionic excitations becomes valid only locally. 
An example for the latter case is a cubic perovskite SrFeOs, "bad" metal which 
is located just slightly below the Mott transition, or CaFe03 having in addition 
a weak charge-density modulation of low-energy holons/doublons. Here, both the 
superexchange interaction accounting for the high-energy charge fluctuations, and 
the low-energy /gapless charge excitations present near the transition, have to be 
considered on equal footing.^ This picture leads to a strong competition between 
the superexchange and double-exchange processes, resulting in orbital disorder, a 
helical spin structure, and small-only Drude weight (quantifying the distance to the 
Mott-transition), as observed in ferratesP^J 

We consider here conventional nearest-neighbor (NN) SE-models, assuming that 
the above criteria 4t 2 /U < t is valid and the local spin-orbital degrees of freedom 
are protected by a charge gap. This is in fact consistent with spinwave measure- 
mentsp*'*™'"^ which can reasonably well be fitted by NN-exchange models in all 
compounds discussed in this section. 

In order to underline a difference between the spin exchange, described by a 
conventional Heisenberg interaction, and that in the orbital sector, we consider first 
the orbital-only models, and move then to the mechanisms which operate in the case 
of full spin-orbital Hamiltonians with different orbital symmetry and spin values. 
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3.1. Orbital-exchange, e g symmetry 

On the cubic lattice, the exchange of e g orbital quantum numbers is described 
by the Hamiltonian 

tfo^frEW' (3-1) 

(v) 

where the pseudospin operators 2t^ = cos 4>^af + sm (fr^af have already been 
used in Ea. ()2-5|) . Here, the orientation of the bond (ij) is specified by the angle 
cp^ = 2yrn/3 with n = 1, 2, 3 for 7 = a, b, c directions, respectively. (For the formal 
derivation of H or i„ consider a spin-polarized version of Ea. l|2-5|) and set Si ■ Sj = 4). 

As pseudospins in (|3-1|) interact antiferromagnetically for all bonds, a staggered 
orbital ordered state is expected to be the ground state of the system. However, 
linear spin-wave theory, when applied to this Hamiltonian, leads to a gapless two- 
dimensional excitation spectrumP*» This results in an apparent instability of the 
ordered state at any finite temperature, an outcome that sounds at least counter- 
intuitive. Actually, the problem is even more severe: by close inspection of the 
orbiton-orbiton interaction corrections, we found that the orbiton self-energy di- 
verges even at zero temperature, manifesting that the linear spin-wave expansion 
about a classical staggered orbital, Neel-like state is not adequate at all. 

The origin of these problems is as follows!^ By symmetry, there are only a 
finite number of directions, one of which will be selected as a principal axis for the 
quadrupole moment. Since this breaks only discrete symmetry, the excitations about 
the ordered state must have a gap. A linear spin wave theory fails however to give the 
gap, because Ea. (|3-l|) acquires a rotational symmetry in the limit of classical orbitals. 
This results in an infinite degeneracy of classical states, and an accidental pseudo- 
Goldstone mode appears, which is however not a symmetry property of the original 
quantum model (|3-1|) . This artificial gapless mode leads to low-energy divergencies 
that arise because the coupling constant for the interaction between orbitons does 
not vanish in the zero momentum limit, as it would happen for a true Goldstone 
mode. Hence the interaction effects are non-perturbative. 

At this point the order-from-disorder mechanism^ comes into play: a particular 
classical state is selected so that the fluctuations about this state maximize the 
quantum energy gain, and a finite gap in the excitation spectra opens, because in 
the ground state of the system the rotational invariance is broken. To explore this 
point explicitly, we calculate quantum corrections to the ground state energy as a 
function of the angle 9 between c-axis and the moment direction. Assuming the 
latter is perpendicular to 6-axis, we rotate globally a pseudospin quantization axes 
as <j z — > a z cos 9 — a x sin 9, a x — > a x cos 6 + a z sin 9, and perform then orbital-wave 
expansion af = 1 — 2a] 04, af ~ a% + a\ around the classical Neel state, where 
the staggered moment is now oriented along the new z direction. As a rotation of 
the quantization axes changes the form of the Hamiltonian, one observes that the 
magnon excitation spectrum has an explicit 0-dependence: 

r 1 1 V2 

Wp(0) = 2A 1 -7 P - — r7 P sin26l - A p (l -cos26l) . (3-2) 
L a/3 J 
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Here, 7 P = (c x + c y )/2, rj p = (c x - c y )/2, A p = (2c z - c x - c y )/6, and c a = cosp a ; 
the energy scale A = 3t 2 /2U. Calculating the zero point magnon energy E{9) = 
— 2 p (-4 — ^u!p(6)), one obtains an effective potential for the staggered moment 
direction (three minima at 9 = <j)^ = 27rn/3), which at small 9 reads as a harmonic 
one: E(9) = const + ^K e jf9 2 , with an effective "spring" constant K e ff = Ak. The 
parameter k is given by 

3 p L (1-7 P ) 1/2 (l-7 P ) d/2J 



The physical meaning of the above calculation is that zero point quantum fluc- 
tuations, generated by interactions between spin waves, are enhanced when the stag- 
gered moment stays about a symmetric position (one of three cubic axes), and this 
leads to the formation of the energy profile of cubic symmetry. A breaking of this dis- 
crete symmetry results then in the magnon gap, which should be about wK^ff/M in 
the harmonic approximation, where an "effective inverse mass" 1/M is of the order of 
the value of the magnon bandwidth W = 2\/2A. More quantitatively, the potential 
E{9) can be associated with an effective uniaxial anisotropy term, \K e ff X^{ij) c °f ° Z ji 
generated by spinwave interactions in the symmetry broken phase. This low energy 
effective anisotropy leads to the gap A = 2^J AK e ff = 2A^[k ~ 0.7^4, stabilizing a 
long-range ordering. The excitation gap compares to the full magnon bandwidth as 
A/W ~ 0.24. Almost the same gap/bandwidth ratio for the model (|3- 1|) was also 
obtained in RefP by using a different method, i.e. the equation of motion. 

The above simple example illustrates a generic feature of orbital-exchange mod- 
els: The interactions on different bonds are competing (like in the three-state Potts 
model), and, very different from the Heisenb erg-like spin interactions, quantum ef- 
fects are of crucial importance even for three-dimensional cubic system. In fact, 
the model (|3-1|) has a finite classical gap and hence no much dynamics in 2D, thus 
fluctuations are more important in 3D. In this way, the orbital degeneracy provides 
a new root to frustrated quantum models in three dimensions, in addition to the 
conventional one driven by geometrical frustration. 

It should be interesting to consider the model (|3-ip for higher dimensional hy- 
percubic lattices, letting the angle be <j>^ = 2irn/d, n = l,...,d. With increasing 
the number of bonds 7 = 1, ...,d the energy potential (having ci-minima as function 
of the moment direction 9) will gradually flatten, and hence the gap will eventually 
close in the limit of infinit dimensions. The ground state and excitations in that 
limit should be very peculiar. 

Finally, by considering a spin-paramagnetic case Si ■ Sj = in Eq. (J2-5p . one 
would arrive again at the orbital model (|3-1|) . which leads to the G-type staggered 
order, different from what is seen in manganites well above T^v- Moreover, spin-bond 
fluctuations in the paramagnetic phase would wash out the three-minima potential, 
hence preventing the orbital order. This indicates again that the orbital transition 
at ~800 K in LaMnOs is primarily not caused by the electronic superexchange. 
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3.2. Orbital-exchange, t2 g symmetry: YTiO^ 

Now, we discuss the i2#-counterpart of the model <|3-lj) . which shows a more 
rich behavior. This is because of the higher, threefold degeneracy, and different 
hopping geometry of a planar-type t2 g orbitals that brings new symmetry elements. 
The orbital order and fluctuations in t^g orbital-exchange model have been studied 
in Refs.^*^ in the context of ferromagnetic YTiC>3. The model reads (in two 
equivalent forms) as: 

H orb = ^- J2 (n ■ rj + \n inj ) (7) (3-4) 

= -=r ^2(n ia n ja + n^n^ + aj/^jaj + ^jaiO^Pj)^ . (3-5) 
1 <*i> 

This is obtained from Ea. (|3-ll|) below by setting S{ ■ Sj = 1/4 (and dropping a 
constant energy shift). The energy E\ = U — 3Jh corresponds to the high-spin 
virtual transition in the spin-polarized state of YTiC>3. In above equations, a ^ f3 
specify the two orbitals active on a given bond direction 7 (see Fig^). For each 
(a/3) -doublet, the Heisenb erg-like pseudospin : ri 7 ^ = [n a —n/})/2, rj 7 ^ = 
= 0^a, and density rS" 1 ^ = n a + rip are introduced. 



C 




Fig. 1. For every bond of the cubic crystal, two out of three t2 g orbitals are equally involved in the 
superexchange and may resonate. The same two orbitals select also a particular component of 
angular momentum. (After RefP'). 

New symmetry elements mentioned above are: (i) a pseudospin interactions are 
of the Heisenberg form; thus, the orbital doublets like to form singlets just like spins 
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do. (m) Apart from an obvious discrete cubic symmetry, the electron density on 
each orbital band is a conserved quantity. Formally, these conservation rules are 
reflected by a possibility of uniform phase transformation of orbiton operators, that 
is, a — > a exp(i(j) a ), which leaves the Hamiltonian invariant. Moreover, as t2g-orbitals 
can hop only along two directions (say, xy-orbital motion is restricted to ab planes), 
the orbital number is conserved on each plane separately. 

The above features make the ground state selection far more complicated than 
in case of e g — orbitals, as it has in fact been indicated long ago. 481 In short (see for 



via the order-from-disorder scenario again. It turns out, however, that in this case 
quantum fluctuations select the body diagonals of the cube as a principal axes for 
the emerging quadrupole order parameter (see FigEJ). The ordered pattern has a 
four-sublattice structure, and the density distribution for the first sublattice (with 
[111] as a principal axis) is described by a function: 



(Similar expressions can easily be obtained for other sublattices by a proper rotation 
of the quantization axes according to Fig|2J). Because of quantum fluctuations, the 
quadrupole moment Q, which controls the degree of orbital elongation, is anoma- 
lously small: Q ~ 0.19 (classically, Q = 1). 



Fig. 2. t2 9 -electron density in the quadrupole ordered state, calculated from Eq. l|M-fi[l . (After 



Surprisingly, not only quadrupole but also a magnetic ordering is equally good 
for the t2 g orbital model. This corresponds to a condensation of complex orbitals 
giving a finite angular momentum, which is again small, mi — 0.19^#- A magnetic 





(3-6) 




RefP). 
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pattern is of similar four-sublattice structure. Further, it turns out that quadrupole 
and magnetic orderings can in fact be continuously transformed to each other - 
using a phase freedom present in the ground state — and a mixed states do appear 
in between. We found that these phase degrees of freedom are the gapless Goldstone 
modes, reflecting the "orbital color" conservation rules discussed above. 

On the technical side, all these features are best captured by a radial gauge 
formalism, applied to the orbital problem in RefP Within this approach, we rep- 
resent the orbiton operators, entering in Ea. (|3-5j) . as a, = ^Pi a e t0ia , thus separating 
a density and phase fluctuations. As a great advantage, this makes it possible to 
single out the amplitude fluctuations (of the short-range order parameters Q, mi), 
responsible for discrete symmetry breaking, from a gapless phase-modes which take 
care of the conservation rules. This way, the ground state condensate wave-function 
was obtained as 

^1,2,3,4(0, <p) = VP^{dyze i{ip+9) ± d zx e^-V ± d xy }. (3-7) 

Here, po <ti 1 determines the amplitude of the local order parameter, while the phases 
ip, 9 fix its physical nature — whether it is of quadrupole or magnetic type. Specifi- 
cally, quadrupole and magnetic orderings are obtained when p = 9 = 0, and <p = it, 
9 = vr/3, respectively. While short-range orbital ordering (a condensate fraction p$) 
is well established at finite temperature via the order-from-disorder mechanism, true 
long-range order (a phase-fixing) sets in at zero temperature only. Slow space-time 
fluctuations of the phases <p, 9 are manifested in a 2D gapless nature of the orbital 
excitations. 

In Eq. (|3-7|) . we recognize the "orbital-angle" picture but, because the order 
parameteres are weak (po = Q/3 ~ 0.06), it represents only a small coherent fraction 
of the wave function; the main spectral weight of fluctuating orbitals is contained 
in a many-body wavefunction that cannot be represented in such a simple classical 
form at all. 

The low-energy behavior of the model is changed once perturbations from lattice 
distortions are included. An important effect is the deviation of the bond angles 
from 180° of an ideal perovskite; this relaxes the orbital-color conservation rules, 
making a weak orbital order possible at finite temperature. Physically, however, 
this temperature is confined to the spin ordering, since the interactions as shown in 
Ea. (|3-5j) are formed only after the spins get fully polarized, while fluctuating spins 
destroy the orbital order which is so fragile even in the ground state. 

A remarkable feature of the SE driven orbital order is that, although cubic sym- 
metry is locally lifted by a small quadrupole moment, the bonds remain perfectly 
the same, as evident from Fig(2J This immediately explains a cubic symmetry of the 
spin-exchange couplings and of the Raman light scattering from orbital fluctua- 
tions,^ — the observations which seem so unusual within a crystal-field picture for 
YTiC>3. This indicates a dominant role of the superexchange mechanism in titanates. 

Orbital excitations in YTiO^. — The superexchange theory predicts the following 
orbital excitation spectrum for YTiOaP^ 

w±(p) = W orb {l - (1 - 2e)(l - 2/) (ti ± k) 2 - 2(e - /)( 7 i ± k)} 1/2 , (3-8) 
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where subscript (±) indicates two orbiton branches, k 2 =72+73? an d 71,2,3 are the 
momentum dependent form-factors 71 (p) = (c x + c y + c z )/3, 72(p) = Vs(c y — c x )/6, 
73(p) = (2c^— c x — c y )/6 with c a = cosp a . Physically, the parameter e ~ 0.2 accounts 
for the many-body corrections stemming from the interactions between orbital waves, 
which stabilize a weak orbital order via the order-from-disorder mechanism. While 
the correction / ~ 0.1 determines the orbital gap: it would be zero within the 
model (|3- 5|) itself, but becomes finite once the orbital-nondiagonal hoppings (induced 
by octahedron tilting in GdFe03-type structure) are included in the calculations. 
Finally, the parameter W or b ~ 2(4t 2 / E±) represents an overall energy scale for the 
orbital fluctuations. By fitting the spin-wave data, 4t 2 /E 1 ~ 60 meV has been 
estimated for YTi03 in Ref.J^ accordingly, W^b ~ 120 meV follows. 

The energy scale 4t 2 /E\ ~ 60 meV is in fact suggested also by a simple estima- 
tion, consider e.g. t ~ 0.2 eV and E\ < 2.5 eV inferred from opticsP^ However, it 
would be a good idea to "measure" this energy from optical experiments, as done 
for manganites."^ 1 At low temperature, when a spin "filtering factor" (Si ■ Sj + 3/4) 
is saturated for the high-spin transition — see first line in Ea. (|3-ll|) — the ground 
state energy (per bond) is: 

4t 2 r,^ . l vM li 4t 2 



E x 



1 1 

Xj H — rum ) 

4 3 



3^i 



(\E \ + 1) , (3-9) 



consisting of a quantum energy Eq ~ —0.214 (per site) calculated in Ref.p^* and a 

constant stemming from a linear terms nj 7 ^ in (|3-llj) . Using now Eq. (|2-8jl . we find 
that Ki ~ 0.81 x (At 2 /Ex), which can directly be determined from the optical carrier 
number N e ff \ once measured. 

Orbital fluctuations in Raman light scattering. — The superexchange energy for 
orbital fluctuations, W or b ~ 120 meV, is apparently consistent with the energy of 
a weak signal in the optical transmission data of Ref.^ Namely, our interpre- 
tation, based on two-orbiton absorption (similar to the two-magnon peak in spin 
systems), gives 2W or b ~ 240 meV for the peak position (without a phonon assist- 
ing the process), as observed. The same characteristic energy follows also for the 
Raman scattering,^ which is derived from a two-orbiton propagator with proper 
matrix elements. However, the most crucial test is provided by a symmetry of the 
orbital states, controlling the polarization dependences. The superexchange theory 
outlined above predicts a broad Raman band (due to a many-body interactions be- 
tween the orbital excitations), with the polarization rules of cubic symmetry. Our 
super exchange-Raman theory is conceptually identical to that for the light scattering 
on spin fluctuations, and, in fact, the observed orbital-Raman lineshapes are very 
similar to those of spin-Raman response in cuprates. On the contrary, we found, 
following the calculations of Ref.J^ that the polarization rules for the lattice driven 
orbital states (|2T0j) in YTi03 strongly disagree with cubic symmetry: The energy 
positions are completely different for the c axis and ab plane polarizations. Such a 
strong anisotropy is imposed by a "broken" symmetry of the lattice: in a crystal- 
field picture, the orbital state is "designed" to fit these distortions (by tuning the 
orbital-angles). 

Comparing the above two models, (|3T|) and ()3-4|) . we see completely different 
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low-energy behavior — while a finite temperature long-range order is protected by a 
gap in former case, no gap is obtained for the t2 g orbitals. Physically, this is related 
to the fact that t2 g triplet accommodates not only the electric quadrupole moment 
but — different from e g doublet — also a true magnetic moment which is a vector 
like spin. It is this dual nature of the tig triplet — a Potts-like quadrupole and 
Heisenberg-like vector — which lies at the origin of rich physics of the model 1)3 • 4 j) . 

3.3. e g spin-orbital model, spin one-half 

Let us move now to the spin-orbital models that describe a simultaneous ex- 
change of both spin and orbital quantum numbers of electrons. We start with the 
SE model for e g holes as in perovskites like KCUF3, neglecting the Hund's rule cor- 
rections for simplicity. On a three-dimensional cubic lattice it takes the formP 

H = J£ $ ■ Sj + - r^\\ - Tj f\ (3-10) 

(ij) 

where J = 4t 2 /U, and t^ 7 ) are the e 9 -pseudospins defined above. The main feature 
of this model — suggested by the very form of Hamiltonian (|3-10|l — is the strong 
interplay between spin and orbital degrees of freedom. It was first recognized in 
Ref.f^ that this simple model contains rather nontrivial physics: the classical Neel 
state is infinitely degenerate in the orbital sector thus frustrating orbital order and 
vice versa; this extra degeneracy must be lifted by some mechanism (identified later- 
on in Ref|23). 

We first notice that the effective spin-exchange constant in this model is definite 
positive for any configuration of orbitals (as r < 1/2), where its value can vary from 
zero to J, depending on the orientation of orbital pseudospins. We therefore expect 
a simple two-sublattice antiferromagnetic, G-type, spin order. There is however a 
problem: a classical G-type ordering has cubic symmetry and can therefore not lift 
the orbital degeneracy. In more formal terms, the spin part (Si ■ Sj + 1/4) of the 
Hamiltonian 1)3-10(1 simply becomes zero in this state for all bonds, so that these 
orbitals effectively do not interact — they are completely uncorrelated. In other 
words, we gain no energy from the orbital interactions. This shows that from the 
point of view of the orbitals the classical Neel state is energetically a very poor. 

The mechanism for developing intersite orbital correlations (and hence to gain 
energy from orbital ordering) must involve a strong deviation in the spin configura- 
tion from the Neel state — a deviation from (Si ■ Sj) = This implies an intrinsic 
tendency of the system to develop low-dimensional spin fluctuations which can most 
effectively be realized by an ordering of elongated 3z 2 — r 2 orbitals [that is, cti = in 
Ea. (|2-l|) ]. In this situation the effective spin interaction is quasi one- dimensional, so 
that spin fluctuations are enhanced as much as possible and large quantum energy is 
gained from the bonds along the 3z 2 — r 2 orbital chains. Since (Si ■ Sj + \) c < 0, the 
effective orbital-exchange constant that follows from 1)3-10)1 is indeed ferromagnetic, 
thus supporting 2>z 2 — r 2 type uniform order. At the same time the cubic symmetry 
is explicitely broken, as fluctuations of spin bonds are different in different directions. 
This leads to a finite splitting of e g — levels, and therefore an orbital gap is generated. 
One can say that in order to stabilize the ground state, a quadrupole order and 
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anisotropic spin fluctuations support and enhance each other — one recognizes here 
the order-from-disorder phenomenon again. 

More quantitatively, the expectation value of the spin-exchange coupling along 
the c-axis, given by strong 3z 2 —r 2 orbital overlap [consider = —1/2 in Ea. (|3-lU|l ]. 
is J c = J, while it is only small in the afr-plane: J a b = J/16. Exchange energy is 
mainly accumulated in c-chains and can be approximated as J c (Si- Sj + \) c + 2J a b(Si- 
Sj + \) a b ~ — 0.16J per site (using (Si ■ Sj) c = 1/4 — In 2 for ID and assuming 
(Si • Sj) a b ~ 0). On the other hand, x 2 — y 2 orbital ordering results in the two- 
dimensional magnetic structure (J a ,b = 9 J/16, J c = 0) with a much smaller energy 
gain ~ — 0.09J. 

From the technical point of view, it is obvious that a conventional expansion 
around the classical Neel state would fail to remove the orbital degeneracy: Only 
quantum spin fluctuations can lead to orbital correlations. This is precisely the rea- 
son why one does not obtain an orbital gap in a linear spin-wave approximation, and 
low-energy singularities appear in the calculations,^ leading to an apparent collapse 
of the spin and orbital order. However, as demonstrated in Ref.p2» the long-range 
orderings are stable against fluctuations, and no "quantum melting" does in fact 
occur. The singularities vanish once quantum spin fluctuations are explicitely taken 
into account in the calculations of the orbiton spectrum. These fluctuations generate 
a finite gap (of the order of J/ 4) for a single orbital as well as for any composite spin- 
orbital excitation. The orbital gap removes the low-energy divergencies, protecting 
the long-range spin order. However, spin order parameter is strongly reduced to 
(S z ) ~ 0.2, due to the quasi-one dimensional structure of spin-orbital correlations. 
Such a spin reduction is generic to spin-orbital models and occurs also in the t2 g 
case, but the mechanism is quite different, as we see below. 

Physically, because of the strong spatial anisotropy of the e g -orbitals, it is im- 
possible to optimize the interactions in all the bonds simultaneously; this results 
in orbital frustration. The frustration is removed here by reducing the effective di- 
mensionality of the interactions, specifying strong and weak bonds in the lattice. 
(We may speak of the "Peierls effect" without the phonons; this is essentially what 
happens in vanadates, too, see later). At the same time tunneling between different 
orbital configurations is suppressed: the spin fluctuations produce an energy gap for 
the rotation of orbitals. A similar mechanism of resolving the frustrations by using 
the orbital degrees of freedom has recently been discussed in Ref.™ for vanadium 
spinels. 

Our last point concerns the temperature scales, T or ,{, and Tat. They are in fact 
both controlled by the same energy, that is the orbital gap A ~ J/4p3 Once the 
quadrupole order is lost at T or b ~ A due to the flat 2D orbital modes, spin order 
will also collapse. Alternatively, thermal destruction of the spin correlations washes 
out the orbital gap. Thus, T or b ~ T/v ~ A in the e g exchange-model alone. [To 
obtain T or b 3> T/v as commonly observed in e g compounds experimentally, the orbital 
frustration should be eliminated by lattice distortions]. In t2 g systems, however, a 
"delay" of T or b, and even T or b <C Tjy, is possible. This is because the t2 g orbitals are 
far more frustrated than the Heisenberg spins are. Such an extreme case is in order 
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to be analyzed now. 

3.4. t2 9 spin-orbital model, spin one-half: LaTiOj, 

We consider first a full structure of the SE Hamiltonian in titanates. Virtual 
charge fluctuation spectra for the Ti-Ti pair is represented by a high-spin transition 
at Ei = U — 3Jh and low-spin ones at energies E% = U — Jh and £3 = U + 2Jh- 
Here, U = A + AB + 3C and Jh = 32? + C are the intraorbital repulsion and 
Hund's coupling in the Ti 2+ excited state, respectively.^ From the optical data 
of Ref.j > one may infer that these transitions are located within the energy range 
from ~ 1 eV to ~ 4 eV. Because of the small spin value, the Hund's splittings 
are apparently less than the linewidth, thus these transitions strongly overlap in 
optics. Indeed, the free-ion value Jh — 0.59 e"V^ for Ti 2+ gives E2 — E\ ~ 1.2 eV, 
compared with a t2 g bandwidth ~ 2 eV. (We believe that Jh is further screened in 
the crystal, just like in manganitea^ ) . Experimentally, the temperature dependence 
of the optical absorption may help to resolve the transition energies, and to fix 
thereby the values of U and Jh in crystal. For YTiOs, we expect that the E\— band 
should increase (at the expense of the low-spin ones) as the sample is cooled down 
developing ferromagnetic correlations. The situation in AF LaTi03 is, however, 
much more delicate because of strong quantum nature of spins in this material (recall 
that the spin-order parameter is anomalously small), and because of the absence of a 
cooperative orbital phase transition. Thus, we expect no sizable thermal effects on the 
spectral weight distribution within the didj— optical multiplet in LaTiC>3. (Optical 
response theory for LaTiC>3, where the quantum effects are of vital importance, is 
still lacking). This is in sharp contrast to manganites, where the classical spin- and 
orbital-orderings lead to a dramatic spectral weight transfer at Tjy and T s t r \—^ 

The above charge fluctuations lead to the SE Hamiltonian^'^ which we repre- 
sent in the following form: 



" = % ^•^ + j)(4j , -i«, W -^ T) ) (3-11) 



+l(a-*-I)(4 , 4". w+ k ) 



\E 3 E 2 )\ l 3 4/3 ij ' 

The spin-exchange constants (which determine the magnon spectra) are given by a 
quantum-mechanical average of the following operator: 



4" = ^ 



\{n + 1-2)4? - \(n - r 3 )B<;f - i(ri - r 2 )(n, + Bj )WJ , (3-12) 

where J = 4t 2 /U. The parameters r n = U/E n take care of the J#-multiplet splitting, 
and r n = 1 in the limit of Jh = 0. One should note that the spin-exchange constant is 
only a fraction of the full energy scale, represented by J, because of the compensation 
between contributions of different charge excitations E n . This is typical when orbital 
degeneracy is present, but more pronounced for t2 g systems where the spin interaction 
may have either sign even in the Jh = limit, see below. 
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The orbital operators A^J , B\- and n\ depend on the bond direction 7, and 
can be represented in terms of constrained particles a, , hi , Cj with rn a + rtj& + n.j c = 1 
corresponding to ti g levels of yz, zx, xy symmetry, respectively. Namely, 

A]f = nianj a + n ib n jb + a\bib^aj + frjajajfy, (3T3) 
•By = nianja + n ib n jf) + a^a^bj + bjaibjaj, 

(c) 

and n\ = ni a +nib, for the pair along the c axis. Similar expressions are obtained for 
the a and b bonds, by replacing (ab)— doublets by (6c) and (ca) pairs, respectively. 
It is also useful to represent A-J and in terms of pseudospins: 

4 7) = 2(n ■ r 3 + Irani) W , B« = 2ft ® r, + i^ry) W , (3-14) 

where 7*; operates on the subspace of the orbital doublet (a, (3)^ active on a 
given 7-bond (as already explained above), while a symbol ® denotes a product 
n ® 75 = 7? r/ + (r+r+ + rr T r)/2. 

At large Jh , the ground state of the Hamiltonian l)3-lip is obviously ferromag- 
netic, imposed by the largest E\— term, and the problem reduces to the model (|3-4j) . 
in which: (i) the orbital wave function is described by Ea. (|3-7|l but, we recall that 
this is only a small condensate fraction; (ii) low-energy excitations are the 2D, gap- 
less, two-branch, Goldstone phase-modes. Concerning the spin excitations, we may 
anticipate some nontrivial things even in a ferromagnetic state. Once a magnon is 
created, it will couple to the orbital phase modes and vice versa. At very large Jh, 
this coupling is most probably of perturbative character but, as Jh is decreased, a 
bound states should form between the spin and orbital Goldstone modes. This is be- 
cause the magnons get softer due to increased contributions of the E% and E% terms. 
Evolution of the excitation spectra, and the nature of quantum phase transition(s) 
with decreasing Jh [at which critical value(s)? of which order?] have not yet been 
addressed so far at all. Needless to say, the finite temperature behavior should be 
also nontrivial because of the 2D modes — this view is supported also by RefP^ 

Looking at the problem from the other endpoint, Jh = 0, E n = U, where the 
ferromagnetic state is certainly lost, one encounters the following Hamiltonian: 

H = 2J^ & ■ S, + \) (tS • r 3 + ^n,-) (7) . (3-15) 

(ii) 

(An unessential energy shift, equal to —J, is not shown here). This model best 
illustrates the complexity of tig orbital physics in perovskites. Its orbital sector, 
even taken alone as in (|3-4jl . is nearly disordered; now, a fluctuating spin bonds will 
introduce strong disorder in the orbital sector, "deadly" affecting the orbital-phase 
modes, and hence the long-range coherence which was already so weak. Vice versa, 
the orbital fluctuations do a similar job in the spin sector; — thus, the bound states 
mentioned above and spin-orbital entanglement are at work in full strength now. 

A while agoP we proposed that, in the ground state, the model (|3-15[) : (i) has 
a weak spin order of G-type which respects a cubic symmetry; (ii) the orbitals are 
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fully disordered. Calculations within the framework of 1/N— expansion, supporting 
this proposal, have been presented in that work. Here, we would like to elaborate 
more on physical ideas that have led to the orbital-liquid picture. 

Obviously, quantum dynamics is crucial to lift a macroscopic degeneracy of clas- 
sical states in the model (|3T5[) . stemming from an infinite number of the "orbital- 
color conservation" rules discussed above. Various classical orbital patterns (like a 
uniform (xy + yz + zx)/y/3, xy orderings, etc.) leave us with Heisenberg spins alone, 
and hence give almost no energy gain and are ruled out. Quasi-ID orbital order like 
in the case of the e g model ()3T0j) is impossible because of a planar geometry of the 
tig orbitals. Yet, the idea of a (dynamical) lowering of the effective dimensionality 
is at work here again, but underlying mechanism is radically different from that in 
Gg case. 

The key point is a possibility to form orbital singlets. Consider, say, the exchange 
pair along the c direction. If both ions are occupied by active orbitals (nf = 

rij = 1), one obtains the interaction of the form 2 J (Si • Sj + 1/4) (rj • fj + 1/4) 
that shows perfect symmetry between spin and orbital pseudospin. The pair has 
sixfold degeneracy in the lowest energy state: both spin-triplet® orbital-singlet and 
spin- singlet® orbital-triplet states gain the same exchange energy equal to — J/2. In 
other words, spin exchange constant may have equally strong ferromagnetic and 
antiferromagnetic nature depending on the symmetry of the orbital wavefunction. 
This violates the classical Goodenough-Kanamori rules, in which ferromagnetic spin 
exchange occurs only at finite Hund's coupling and hence is smaller by factor of 
Jh/U. In this respect, ti g superexchange clearly differs from the e g model (|3T0|) . in 
which the spin-exchange interaction is positively definite because no orbital singlets 
can be formed in that case. 

When such tig pairs form ID chain, one obtains a model which has been in- 
vestigated under the name SU(4) modelP^'^ A large amount of quantum energy 
(— 0.41J per site) is gained in this model due to resonance between the local configu- 
rations spin-triplet® orbital-singlet and spin-singlet® orbital-triplet. As a result of this 
resonance, low-energy excitations are of composite spin-orbital nature. In a cubic 
lattice, the situation is more complicated, as SU(4) spin-orbital resonance along one 
direction necessarily frustrates interactions in the remaining two directions which 
require different orbital pairs (see Fig^). Given that SU(4) chain physics so ideally 
captures the correlations, one is nevertheless tempted to consider a "trial" state: the 
xy orbital is empty, while xz/yz doublets (together with spins) form SU(4) chains 
along the c axis — a kind of spin-orbital nematic, with a pronounced directionality 
of the correlations. Accounting for the energy- lost on a "discriminated" (classical) 
ab plane bonds on a mean-field level (J/8 per site), we obtain Eq = —0.29 J for our 
trial state, which is by far better than any static orbital state, and also better than 
the ferromagnetic state with fluctuating orbitals (Eq = — 0.214 Once the xy 
orbital is suppressed in our trial state, the interchain couplings read as 

H\f ] = - S 3 + \)(\±t*)(\±tZ), (3-16) 

where ± sign refers to a/b bond directions. In the ground state, these couplings 
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may induce a weak ordering (staggered between the chains) in both sectors, which, 
however, should not affect much the intrachain SU(4) physics, by analogy with ~1D 
spin systems.^ (This should be an interesting point to consider). 

An assumption made is that a quadrupole order parameter, Q = n a + n& — 2n c , 
responsible for the chain structure, is stabilized by order-from-disorder as in case of e g 
quadrupoles in (j3-10j) . or as it happens in the spin-one model for vanadates (see later). 
However, the t2 g quadrupole is highly quantum object, as wee have seen above in 
the context of YTiC>3 , and it is hard to imagine that the above structure will survive 
against the xy orbital intervention, that is, "cutting" the SU(4 ) chains in small pieces 
and their orientational disorder. One may therefore think of a liquid of ,, SU(4)- 
quadruplets" (a minimal building block to form a spin-orbital sing 

lePJEJ). This 

way, one arrives at an intuitive picture of dynamical patterns where the local physics 
is governed by short-range SU(4) correlations, like in quantum dimer models. As a 
crude attempt to capture the local SU(4) physics, we applied the 1/N— expansion 
to the model (|3-15|) . introducing a bond-operator of mixed spin-orbital nature. A 
quadrupole disordered state was found indeed lower in energy (Eq ~ -0.33 J) 71 than 
a nematic state just discussed. As a 1/N— expansion usually underestimates the 
correlations, we think that a quadrupole disordered state best optimizes the overall 
quantum energy of the t2 g spin-orbital superexchange. An additional energy profit 
is due to the involvement of all three orbital flavors — a natural way of improving 
a quantum energy gain. The nature of orbital excitations is the most fundamental 
problem. Tentatively, we believe that a pseudogap must be present protecting a 
liquid state; this has already been indicated in RefP (as an orbital gap, stemming 
formally from the pairing effects within 1/N— expansion). 

An important point is that spins and orbitals in the 3D model (|3-15|) are not 
equivalent. In the spin sector, the Heisenberg interactions on different bonds coop- 
erate supporting the spin long-range order (albeit very weak) in the ground state. 
It is the orbital frustration which brings about an unusual quantum physics in a 3D 
system. When orbitals are disordered, the expectation value of the spin-exchange 
constant, given by Eq. ()3-12|) . is of AF sign at small Jh, supporting a weak spin-G 
order in the ground state, on top of underlying quantum spin-orbital fluctuations. 
Important to note is that the local SU(4 ) physics is well compatible with a weak spin 
staggering. The main ingredient of the theory of RefP-^ is a local SU(4) resonance, 
which operates on the scale of J and lifts the orbital degeneracy without symmetry 
breaking. A remote analogy can be drawn with a dynamical JT physics: — the role 
of phonons are played here by spin fluctuations, and an entangled SU(4) spin-orbital 
motion is a kind of vibronic state but living on the bonds. While orbital-lattice 
vibronic states are suppressed by classical structural transitions, the orbital-spin 
SU(4) resonance may survive in a lattice due to quantum nature of spins one-half 
and orbital frustration, and may lead to the orbital disorder in the 3D lattice — this 
is the underlying idea. 

A weak staggering of spins (while the orbitals are disordered) is due to the 
spin-orbital asymmetry for the 3D lattice. The Hund's coupling Jh brings about 
yet another asymmetry between the two sectors, but this is now in favor of spin 
ferromagnetism. Jh changes a balance between two different (AF and F) spin corre- 
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lations within a SU(4) resonance, and, eventually, a ferromagnetic state with a weak 
3D quadrupole order (Figj^J is stabilized. Are there any other phases in between? 
Our tentative answer is "yes", and the best candidate is the spin-orbital nematic 
discussed above. This state enjoys a fully developed SU(4) physics along the c di- 
rection, supported by orbital quadrupole ordering (xy orbital-selection). The xy 
orbital gap, induced in such a state by Jh in collaboration with order- from-disorder 
effect, is still to be quantified theoretically. In this intermediate phase, spin and 
xz/yz doublet correlations are both AF within the planes (see Eci l3T6[) . but differ- 
ent along the SU(4) chains: more ferro (than AF) for spins, and the other way round 
in the orbital sector. Thus, we predict an intermediate phase with a weak spin-C 
and orbital-G order parameters. Our overall picture is that of the three competing 
phases: (I) spin-ferro and orbitals as in Fig|2 stable at large Jh', (II) spin-C, dou- 
blets xz/yz are staggered, xy occupation is less than 1/3; (III) spin-G/orbital- liquid 
at small Jh- From the experience in vanadates (see next section), we suspect that a 
tight competition between these states may occur for realistic Jh values. The first 
(last) states are the candidates for YTi03 (LaTiOs); it should be a good idea looking 
for the intermediate one at compositions or compounds "in between". Needless to 
say, that all these three states are highly anomalous (compare with 3D Heisenberg 
systems), because the classical orderings here are just a secondary effects on top of 
underlying SU(4) quantum fluctuations (or of pure orbital ones at large Jh)- 

Physically, Jh— tuning is difficult but can be somewhat mimicked by a variation 
of the Ti-O-Ti bond angle 6 (e.g., by pressure). A deviation of it from 180% gives 
an additional term in the spin-exchange through the small t^g — e g overlap as it has 
been pointed out in RefP'' According to Ref.J™ this term supports ferromagnetism 
equally in all three directions (different from RefP^). Thus, such a term: —J' Si ■ Sj 
with J' oc sin 2 does not break a cubic symmetry itself, and hence may perfectly 
drive the above phase transitions. A pronounced quantum nature of the competing 
phases (because of quantum orbitals) may lead to a "soft" character of transitions, as 
suggested in RefP Yet another explanation, based on classical orbital description, 
has been proposed in Ref.f^ predicting the spin- A phase as an intermediate state. 
Thus, the predictions of a quantum and classical orbital pictures are very different: 
the spin-C versus the spin- A type intermediate state, respectively. This offers a nice 
opportunity to discriminate between the electronic and lattice mechanisms of lifting 
the orbital degeneracy in titanates. 

Summarizing, the Hamiltonians l]3-llj) and (|3-15|) are the big " puzzles" , provid- 
ing a very interesting playground for theory. In particular, the phase transitions 
driven by Jh are very intriguing. Concerning again the relation to the titanites: 
While the most delicate and interesting low-energy problems are (unfortunately) 
eliminated by weak perturbations like lattice distortions, the major physics — a 
local SU(4) resonance — should be still intact in LaTi03. This view provides a 
hitherto unique explanation for: (i) an anomalous spin reduction (due to a quan- 
tum magnons involved in the spin-orbital resonance, see RefP); (it) the absence of 
a cooperative structural transition (the orbital liquid has no degeneracy, hence no 
JT instability at small coupling); (iii) nearly ideal cubic symmetry of the spin and 
Raman responses in both LaTi03 and YTi03 (in full accord with our theory). The 
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identification of the predicted intermediate spin-C phase is a challenge for future 
experiments. 

3.5. tig spin-orbital model, spin one: LaVO^ 

In the model for titanites, a quantum nature of spins one-half was essential; to 
make this point more explicit, we consider now similar model but with the higher 
spin, 5 = 1. Apart from its direct relevance to pseudocubic vanadates AVO3, the 
model provides yet another interesting mechanism of lifting the degeneracy by SE 
interactions: here, the formation of the quantum orbital chains is the best solution!^ 

The interactions between 5 = 1 spins of V 3+ ions arise from the virtual exci- 



tations d 2 d 2 - 



djdp and the hopping t is allowed only between two out of three 



tig orbitals, just as in titanites. The df excited state may be either (1) a high-spin 
A2 state, or one of a low-spin ones: (ii) the degenerate 2 E and 2 T\ states, or (Hi) 
a 2 T2 level. The excitation energies are E\ = U — 3Jh, E2 = U and E3 = U + 2 Jh, 
respectively,^ where U = A + AB + 3C and Jh = 3B + C. For the free ion V 2+ , one 
has Jh — 0.64 e\fl^ but this should be screened in crystal to ~ 0.5 eV as suggested 
m RefP Correspondingly, the SE Hamiltonian consists of three contributions, like 
in Eq. ()3-llj) . but a different form as obtained in RefP is more instructive here: 
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1 ' Bjj and introduced previously in Eas. (|3-13j) - ()3-14j) . 
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Here J = At 2 /U, as usual. The coefficients R = U jE\ = 1/(1 — 3n) and r 
1/(1 + 27/) with n = Jh /U take care of the Jh — multiplet structure. 

If we neglect the Hund's splitting of the excited states (consider r\ - 
the Hamiltonian (|3-17j) reduces to: 



H = jY,\(S t -S 3 + l){n-r j + l -n l n j ) 



(7) 



(3-18) 
(3-19) 

: U/E 3 = 

limit), 
(3-20) 



where a constant energy of — 2J per site is neglected. This result should be com- 
pared with corresponding limit in the d 1 case, Ea. (|3-15|) . One observes different spin 
structures: \(§i • Sj + 1) is obtained for vanadium ions instead of 2(5, • Sj + i) for 
spins one-half of Ti 3+ . The difference in spin values can in fact be accounted for in 
general form as (5j ■ Sj + 5 2 )/25 2 . It is important also to note that we have two 
electrons per V 3+ ion; one therefore has a different constraint equation for orbiton 
densities n, a + Uj& + rij c = 2. 
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It is instructive to start again with a single bond along the c-axis. A crucial 
observation is that the lowest energy of —J/2 is obtained when the spins are fer- 
romagnetic, and the orbitals a and b form a singlet, with (fi ■ Tj)^ = — §. Spin 
singlet® orbital triplet level is higher (at —J/4). This is in sharp contrast to the 
S = 1/2 case, where the spin singlet® orbital triplet and the spin triplet® orbital sin- 
glet configurations are degenerate, resulting in a strong quantum resonance between 
them as it happens in titanates. Thus, ferromagnetic interactions are favored due to 
a local orbital singlet made of a and b orbitals. Dominance of high spin configuration 
reflects simply the fact that the spin part of the interaction, that is (5j • Sj + S 2 )/2S 2 , 
is equal to 1 for a ferromagnetic configuration, while vanishing in spin-singlet sector 
(as —1/S) in the limit of large spins. In order to form ab— orbital singlet on the 
bond along c axis, the condition nf* = = 1 must be fulfilled (no pseudospin 
can be formed otherwise). This implies that the second electron on both sites has 
to go to an inactive (that is xy) orbital. Thus we arrive at the following picture for 
the superexchange bond in c direction: (i) spins are aligned ferromagnetically, (ii) 
one electron on each site occupies either a or b orbital forming a SU(2) invariant 
orbital pseudospins that bind into the orbital singlet, (Hi) the xy orbital obtains a 
stabilization energy of about — J/2 (the energy required to break ab— orbital singlet) 
and accommodates a remaining, second electron. 

Formation of one- dimensional orbital chains. — If the high spin state of the given 
pair is so stable, why does then a whole crystal not form uniform ferromagnet? That 
is because each site has two electrons, and an orbital that is inactive in particular 
(ferromagnetic bond) direction, induces in fact an antiferromagnetic coupling in the 
other two directions. Thus spin interactions are strongly ferromagnetic (supported 
by orbital singlets) in one direction, while the other bonds are antiferromagnetic. As 
all directions are a priori equivalent in a cubic lattice, we again arrive at the problem 
of "orbital frustration" common to all spin-orbital models on high-symmetry lattices. 
The solution of this problem here is as follows. As the spin-orbital resonance like 
in titanates is suppressed in the present case of large spin 5=1, quantum energy 
can be gained mainly from the orbital sector. This implies that a particular classical 
spin configuration may be picked up which maximizes the energy gain from orbital 
fluctuations. Indeed, orbital singlets (with ni a + = 1) may form on the bonds 
parallel to the c-axis, exploiting fully the SU(2) orbital physics in one direction, while 
the second electron occupies the third ti g orbital (rij c = 1), controlling spin interac- 
tions in the aft-planes. Thus one arrives at spin order of the C-type [ferromagnetic 
chains along c-axis which stagger within afr-planes], which best comprises a mixture 
of ferromagnetic (driven by the orbital singlets) and antiferromagnetic (induced by 
the electron residing on the static orbital) interactions. This is an analogy of the 
intermediate phase that we introduced for titanites above; here, it is much more 
stable because of the large spin value. 

Once the C-type spin structure and simultaneous selection among the orbitals 
(fluctuating a, b orbitals, and stable c orbital located at lower energy) is made by 
spontaneous breaking of the cubic symmetry, the superexchange Hamiltonian can 
be simplified. We may set now n,- ic = l,nj a + = 1, and introduce pseudospin f 
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which operates within the (a, b) orbital doublet exclusively. 

We focus first on orbital sector as quantum dynamics of the system is dominated 
by the orbital pseudospins r = | rather than by large spins 5 = 1. In the classical 
C-type spin state, (Si ■ Sj + 1) is equal 2 along the c-direction while it is zero on 
a6-plane bonds. In this approximation, orbital interactions in the model l)3T7JI are 
given by (2J-^ +K^') on c-axis bonds, while on a6-plane bonds only the term 
contributes which is small. Expressing A\f and B\f operators in Ens. (l3^ - (l3~T31) 
via pseudospins r, one arrives at the following orbital HamiltonianP^ 

H orb = J orb ^(T l -r,) + J^ rb Yl T i T h ( 3 " 21 ) 

where J orb = JR and J^ rb = Jr](R + r)/2. As their ratio is small, J^/Jorb = 
7/(1 — 5r]r/2) (that is about only 0.1 for the realistic values of parameter 77 = Jh/U 
for vanadates), we obtain one-dimensional orbital pseudospin chains coupled only 
weakly to each other. Orbital excitations in the model ()3-21j) are mostly propagating 
along c-chain directions. Their spectrum can be calculated, e.g., within a linear spin- 
wave approximation, assuming a weak orbital order due to interchain coupling 
One indeed obtains the one-dimensional orbital-wave spectrum; 39 ' 



lj p = ^A2 + J* rb sm 2 p z , (3-22) 

which shows the gap A = J{n(R + r)[2R + r/(R + r)]} 1 / 2 at p z = ir. The orbital gap 
A is small and grows with increasing Hund's coupling as oc J ^fr\. 

Alternatively, one can use the Jordan- Wigner fermionic representation to de- 
scribe quasi one-dimensional orbital dynamics, as suggested in RefP^ One obtains 
the ID orbiton dispersion: 



£k = y/h 2 + Jl b cos 2 k z , (3-23) 

where h = ^rJ^ rh is the ordering field stemming from interchain interactions. The 
staggered orbital order parameter r = |(t*)|, determined self-consistently from r = 
J2k(h/2 £ k) tanh(£fc/2T), is small, and orbitals strongly fluctuate even at T = 0. 

The underlying ID-orbital dynamics have an important consequences on spin 
interactions which control spinwave dispersions. In the spin sector, we obtain inter- 
actions J c (Si ■ Sj) and J a b(Si ■ Sj), with the spin-exchange constants following from 
Ea. (|3T8|) . The result is given by orbital pseudospin correlations: 



j c = j - 

2 

Jab = 1 



;i + 2r/i?)<f i • Tj + i) - W {n ® Tj + i> - r/fl] , (3-24) 
1 - ri(i? + r) + (1 + 2r]R - 7/r)(fi ® Tj + -) 



While in-plane antiferromagnetic couplings are mostly determined by the classical 
contribution of xy orbitals (first term in J ab ), the exchange constant along the 
axis has substantial quantum contribution represented by the first term in J c . This 
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contribution is of negative sign due to the orbital singlet correlations along chains. 
The pseudospin expectation values in ()3-24j) can be estimated by using either the 
Jordan- Wigner fermionic-orbiton representation^ or within orbital-wave approxi- 
mation. 39 ' In both cases, one observes that the ferromagnetic coupling along the 
c-axis is strongly enhanced by orbital fluctuations. For a realistic values of r] ~ 0.12, 
one obtains — J c ~ J a b ~ J/5. This result is qualitatively different from that expected 
from the Goodenough-Kanamori rules. Indeed, in that classical picture with fully 
ordered orbitals one would obtain instead the smaller value — J c ~ 2r]RJ a b ~ Jab/2- 
Comparing all the coupling constants in both spin and orbital sectors, one ob- 
serves that Heisenberg-like orbital dynamics has the largest energy scale, J orb — JR^ 
thus dominating the physics of the present spin-orbital model. The overall picture is 
that "cubic frustration" is resolved here by the formation of the orbital chains with 
Heisenberg dynamics, from which a large quantum energy is gained. This is sim- 
ilar to the order-from-disorder scenario for the e g spin-orbital model (|3- 1U() . where 
classical 3z 2 — r 2 orbital order results in quasi one-dimensional spin chains. In the 
present t<i g orbital model with large classical spins, the role of spins and orbitals are 
just interchanged. 

As argued in Ref.pS) the above scenario may explain the C-AF type spin order 
in LaV03. 60 ' A structural transition that follows the onset of magnetic order in 
this compound is also natural in the superexchange model. Indeed, C-type spin 
ordering and formation of the pseudospin orbital chains are intimately connected 
and support each other. Thus, the ordering process in the orbital sector — the 
stabilization of xy orbital which absorbs one electron, and lifting the remaining 
degeneracy of xz / yz-doublet via quantum fluctuations — occurs cooperatively with 
the evolution of the C-type magnetic order. In more formal terms, the classical order 
parameters Q or b = (2n c — n a — rib) (^y-orbital selection) and Q sp = ((Si ■ Sj) c — (Si • 
Sj) a b)/2 (spin-bond selection) act cooperatively to break the cubic symmetry, while 
the energetics is driven by a quantum energy released by xz/yz orbital fluctuations. 
Obviously, this picture would break down if the JT-coupling dominates — nothing 
but a conventional 3D-classical ordering of xz/y z-doublets at some T s t r , independent 
of spin correlations, would take place in that case. 

A pronounced anisotropy and temperature dependence of the optical conductiv- 
ity observed in LaVOj^J also find a quantitative description 15 ' within this theory, in 
which quantum orbital fluctuations play the key role. Interestingly, the JT picture 
has also been applied to the problem, see RefP^ Based on first-principle calcu- 
lations, a JT binding energy Ejt ~ 27 meV has been obtained. Consequently, a 
large orbital splitting (= 4Ejt) suppresses the quantum nature of orbitals and, as 
a result, one obtains that the optical spectral weghts along the c and a, b directions 
are almost the same, I c /I a b — 1) contradicting the experiment which shows strong 
polarization dependence. In our view, the optical experiments^ clearly indicate 
that JT coupling is much smaller than estimated in Ref.p^ and are thus not suf- 
ficient to lock-in the orbitals in LaV03. At this point, we are faced again with a 
problem of the principal importance: Why do the first-principle calculations always 
predict a large orbital splittings? should those numbers be directly used as an input 
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parameters in the many-body Hamiltonians, suppressing thereby all the interesting 
orbital physics? These questions remain puzzling. 

Based on the above theory, and by analogy with titanites, we expect that the Ra- 
man light scattering on orbital fluctuations in vanadates should be visible as a broad 
band in the energy range of two-orbiton excitations, 2J or b, of the quantum orbital 
chains. Given a single-orbiton energy scale J or b = JR ~ 60 — 70 meV (which follows 
from the fit of opticaP^ and neutron scattering^ data) , we predict a broad Raman 
band centered near ~ 120 meV. This energy is smaller than in titanites because 
of the low-dimensionality of orbital dynamics, and falls in the range of two-phonon 
spectra. However, the orbital-Raman band is expected to have very large width (a 
large fraction of J rb), as orbitons strongly scatter on spin fluctuations (as well on 
phonons, of course). More specifically, we observe that both thermal and quantum 
fluctuations of the spin-bond operator, (Si ■ Sj + 1) in Ea. (|3-17|) or (|3-20f> . introduce 
an effective disorder in the orbital sector. In other words, the orbital Hamiltonian 
(|3-21j) obtains a strong bond-disorder in its coupling constants, hence the expected 
incoherence and broadness of the orbital Raman response. This feature, as well as a 
specific temperature and polarization dependences should be helpful to discriminate 
two-orbiton band from two-phonon (typically, sharply structured) response. On the 
theory side, the light scattering problem in the full spin-orbital model ()3-17j) need to 
be solved in order to figure out the lineshape, temperature dependence, etc. Alter- 
natively, a crystal-field and/or first-principle calculations would be helpful, in order 
to test further the "lattice-distortion" picture for vanadates, — in that case, the 
Raman-band frequencies would be dictated by the on-site level structure. 

Let us conclude this section, devoted to the spin-orbital models: Orbital frus- 
tration is a generic feature of the superexchange interactions in oxides, and leads 
to a large manifold of competing states. Reduction of the effective dimensionality 
by the formation of quantum spin- or orbital-chains (depending on which sector is 
"more quantum") is one way of resolving the frustrations. In e g spin-orbital mod- 
els, the order-from-disorder mechanism completes the job, generating a finite orbital 
gap below which a classical description becomes valid. In superexchange models for 
titanites, where both spins and orbitals are of quantum nature, a composite spin- 
orbital bound states may develop in the ground state, lifting the orbital degeneracy 
without breaking the cubic symmetry. The nature of such a quantum orbital-liquid 
and its elementary excitations is of great theoretical interest, regardless to which 
extend it is "spoiled" in real materials. The interplay between the SE interactions, 
dynamical JT coupling, and extrinsic lattice distortions is a further step towards a 
quantitative description of the electronic phase selection in titanites. 

§4. Competition between superexchange and lattice effects: YVO3 

As an example of the phase control by a competing superexchange, lattice, tem- 
perature and doping effects, we consider now a slightly modified version of the S = 1 
model (|3-17|) for vanadates, by adding there the following Ising-like term which op- 
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erates on the bonds along c direction!^ 

H v = ~ v E T * T r ^ 

fo')ll c 

This describes a ferromagnetic (V > 0) orbital coupling, which competes with the 
orbital-Heisenberg J or b— term in the Hamiltonian (|3-21|) . Physically, this term stands 
for the effect of GdFeC>3 type distortions including A-0 covalency, which prefer a 
stacking ("ferromagnetic" alignment) of orbitals along c direction.^ This effect 
gradually increases as the size of A-site ion decreases and may therefore be partic- 
ularly important for YVO3. Interesting to note that the V— term makes the orbital 
interactions of xyz-type with z-coupling being smaller, hence driving the orbital sec- 
tor towards a more disordered xy model limit (motivating the use of Jordan- Wigner 
representation for orbitals 15 '). However, when increased above a critical strength, 
the V— term favors a classically ordered orbitals. 

Three competing phases can be identified for the modified, (|3-17|) plus (|4-1|) . 
model, (i) The first one, which is stable at moderate V < J, was that just considered 
above: C-type spin order with fluctuating orbitals (on top of a very weak staggered 
order) as in LaVC>3. This state is driven by a cooperative action of orbital-singlet 
correlations and J#- terms, (ii) At large V > J, the orbitals order ferromagnetically 
along the c-axis, enforcing the G-type spin order, as observed in YVO3. (Hi) Yet 
there is a third state, which is "hidden" in the SE-model and may become a ground 
state when both Jjj and V are below some critical values. This is a valence bond solid 
(VBS) state, which is derived by a dimerization of the Heisenberg orbital chains, by 
using the spin-bond operators (Si ■ Sj + 1) as an "external" dimerization field. (We 
already mentioned this operator as a "disorder field" for the orbital excitations). 

Specifically, consider the limit of r\ = 0, V = 0, where we are left with the 
model ()3-2(j[) alone. It is easy to see that we can gain more orbital quantum energy 
by choosing the following spin structure: On every second c-axis bond, spins are 
ferromagnetic, (Sj • Sj + 1) = 2, while they are antiparallel in all other bonds giving 
(Si ■ Sj + 1) ~ 0. Consequently, the orbital singlets are fully developed on ferro-spin 
bonds, gaining quantum energy —J/2 and supporting the high-spin state assumed. 
On the other hand, the orbitals are completely decoupled on AF-spin bonds. As 
a result, the expectation value (f| • fj) vanishes on these "weak" bonds, thus the 
spin exchange constant Jf eak = J/8 (consider Eq. (|3-24l) for uncorrelated orbitals 
and the r/ = limit) turns out to be positive, consistent with antiferromagnetic spin 
alignment on these bonds. Such a self-organized spin-orbital dimer state is lower in 
energy [—J /A per site] than a uniform Heisenberg orbital chains with C-type spin 
structure [(1/2 - In 2) J ~ -0.19J]. Thus, the VBS state (building blocks are the 
decoupled orbital dimers with total spin 2) is a competing state at small rj values, 
as has been noticed by several &\\th.ors^^*^Il^i®^ 

4.1. Phase diagram 

Interplay between the above three spin-orbital structures has been investigated 
numerically in RefP within the DMRG method. The thermodynamic properties 
and interesting evolution of correlation functions have also been studied by using 
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the finite-temperature version of DMRG.BDBB Though these methods are restricted 
to the one-dimensional version of the model, they give rigorous results and are in 
fact quite well justified, because the essential physics is governed by strong dynamics 
within the ~1D orbital chains, while a weak interchain couplings can be implemented 
on a classical mean-field level. 11 : 




Fig. 3. Phase diagram in the rj-V plane. The VBS state consists of orbital-singlet dimers with 
spin 2, and the spins of different dimers are weakly coupled antiferromagnetically. In the phase 
C, a uniform orbital-chain is restored, while the spins are aligned ferromagnetically along the 
chain. The phase G is the spin- AF/orbital- ferromagnetic state, stabilized by the Ising interaction 
between the orbitals originating from the GdFeOs-type distortions. All the phase transitions 
are of first-order. (After RefP*). 



The ground state phase diagram in the t)-V plane, obtained from the DMRG 
study^ is shown in FigOl There are three distinct phases in this figure. For 
small n and V, the VBS state is stabilized, which is driven either to the orbital- 
ferromagnetic/spin- AF phase (called G) with the increase of V, or to the orbital- 
AF/spin- ferromagnetic one (called C) with the increase of the Hund's coupling rj. 
The critical value r] c (V = 0) ~ 0.11 for the latter transition perfectly agrees with 
the earlier result inferred from the finite-temperature DMRG-studyj^ and is just 
slightly below the realistic values in vanadates.^' This indicates the proximity of 
VBS state in vanadates, on which we elaborate later-on. Note that when the Hund's 
coupling is slightly above the critical value, the ground state is spin-ferromagnetic for 
small V, but for intermediate values of V the VBS-state is stabilized. Stabilization 
of the orbital disordered state by finite V interaction (induced by lattice distortion) 
is a remarkable result. The physics behind this observation is that — as already 
mentioned — the interaction (|4-1|) introduces a frustration into the orbital sector, 
competing with antiferro-type alignment of orbitals in the SE model. 

The V— driven phase transition from the spin-C to spin-G ground state, obtained 
in these calculations describes physically why an increased GdFeC>3-type distortion 
promotes a spin-staggering along the c direction in YVO3, while less distorted LaV03 
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has the spin-C structure. This study also suggests that V ~ J in vanadates. As 
J ~ 40 meV in these compounds we suspect that the effect of GdFeC>3-type 
distortions on t2 g orbitals is roughly of this scale in perovskites in general, includ- 
ing titanites. Being comparable with the SE scale J or b ~ 60 meV in vanadates, 
the V— term is important to stabilize the spin-G phase, but it is not sufficient to 
lock-in the orbitals in titanates, where the three dimensional exchange fluctuations 
are characterized by larger energies of the order of W or b ~ 120 meV (see previous 
Section). 

4.2. Entropy driven spin-orbital dimerization 

YVO3 is a remarkable material, in the sense that its spin-G type ground state 
is very fragile and readily changes to the spin-C state via the phase transition, 
driven either by temperature (at 77 K) or by small doping. 66 ' This finds a natural 
explanation within the present theory in terms of competing phases, assuming that 
YVO3 is located in the phase G near the tricritical point (see the phase diagram in 
Fig|3J), and thus close to the spin-C and VBS states. 

This view is strongly supported by the neutron scattering data of RefP" This 
experiment revealed several anomalies, indicating that the spin-C phase above 77 K 
is itself very unusual: (i) substantial modulations of the spin couplings along the c 
direction, (ii) ferromagnetic interactions are stronger than oh— plane AF ones, (iii) 
anomalously small ordered moment, which (iv) selects the ab— plane (different from 
the easy c axis found just below 77 K or in the C-phase of LaVO^). All these features 
have coherently been explained in RefP^» in terms of underlying quantum dynamics 
of the spin-orbital chains and their entropy- driven dimerizationP^ Physics behind 
the last point is as follows. 

Because the Hund's coupling parameter is close to the critical one, there is strong 
competition between uniform (C) and dimerized (VBS) states, and this may affect 
thermodynamic properties for entropy reasons. The point is that the dimer state 
contains "weak" spin bonds: spin interaction between different dimers is small when 
77 is near the critical value (when the dimerization amplitude is not saturated, weak 
bonds are ferromagnetic but much weaker than strong ferromagnetic bonds within 
the dimers). Therefore, the spin entropy of an individual dimer with total spin 2, 
that is In 5, is released. The gain of spin entropy due to the dimerization lowers 
the free energy F = (H) — TS and may stabilize a dimerized state with alternating 
weak and strong ferro-bonds along c-axis. In other words, the dimerization of spin- 
orbital chains occurs due to the orbital Peierls effect, in which thermal fluctuations 
of the spin bond-operator (Si ■ Sj) c play the role of lattice degrees of freedom, while 
the critical behavior of the Heisenb erg-like orbital chains is a driving force. As the 
dimerization is of electronic origin, and is not as complete as in the VBS state itself, 
concomitant lattice distortions are expected to be small. 

These qualitative arguments have been substantiated by numerical studies using 
the finite-temperature DMRG method.^^^ The explicit calculations of the entropy, 
evolution of the dimer correlations upon increasing temperature, and anomalous 
behavior of spin and orbital susceptibilities can be found in these papers. 
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4.3. Doping control 

The energy difference between G and G type spin orderings in YVO3 is very 
small, E c -E G ~ 0.1J ~ 4 meV only® J® 

for realistic values of 77 and V. There- 
fore, the G-type ground state of YVO3 can easily change to the C-type upon small 
perturbations. For instance, pressure may reduce the GdFeC>3 distortions hence the 
V— interaction, triggering the first-order phase transition described above. Injec- 
tion of the charge carries is the another possibility, which we discuss now (see also 
Ref.®). 

We need to compare a kinetic energy gains in above-mentioned competing phases. 
It is easy to see that the charge-carriers strongly favor spin-G phase, as they can 
freely move along the ferromagnetic spin chains in the c direction of the G phase, 
while a hole-motion is frustrated in the spin-G phase in all directions. In more detail, 
the afr-plane motion is degraded in both phases equally by a classical AF order via 
the famous "double-exchange factor" cos | with 8 = tt for antiparallel spins. Thus, 
let us focus on the c direction. In the G-phase, a hole-motion is disfavored again 
by antiparallel spins (9 C = tt), but the spin-G phase with 9 C = is not "resistive" 
at all. Now, let us look at the orbital sector. The orbitals are fully aligned in the 
spin-G phase and thus introduce no frustration. Our crucial observation is that the 
orbitals are not resistive to the hole motion in the spin-G phase, too. The point 
is that a doped-hole in the spin-G phase can be regarded as a fermionic-holon of 
the quasi-lD orbital chains, and its motion is not frustrated by orbitals at all. This 
is because of the orbital-charge separation just like in case of a hole motion in ID 
Heisenberg chains!® As a result, the c direction is fully transparent in the spin- 
G phase (both in spin and orbital sectors), and doped holes gain a kinetic energy 
K = 2tx per site at small x. (On contrast, a hole motion is strongly disfavored in 
spin-G phase because of AF-alignment of large spins S = 1, as mentioned). Even for 
the doping level as small as x = 0.02, this gives an energy gain of about 8 meV (at 
t ~ 0.2 eV) for the G phase, thus making it lower in energy than the G-phase. This 
is presicely what is observed.® The underlying quantum orbital chains in the spin- 
G phase® are of crucial importance here: a static configuration of the staggered 
xz/yz orbitals would discourage a hole-motion. In other words, fluctuating orbitals 
not only support the ferromagnetic spins, but also well accommodate doped holes 
because of orbital-charge separation. The quantum orbitals and a doping induced 
double-exchange act cooperatively to pick up the spin-G phase as a ground state. 

The "holon-like" behavior of the doped- holes implies a quasi one-dimensional 
charge transport in doped vanadates, namely, a much larger and strongly tempera- 
ture dependent low-energy optical intensity (within the Mott-gap) along the c axis, 
in comparison to that in the afc-plane polarization. The situation would be quite 
different in case of classical JT orbitals: A staggered 3D-pattern of static orbitals 
frustrates hole motion in all three directions, and only a moderate anisotropy of the 
low-energy spectral weights is expected. This way, optical experiments in a doped 
compounds have a potential to discriminate between the classical JT picture and 
quantum orbitals, just like in the case of pure LaV03 as discussed above. 

The doping-induced transition from the G to the G phase must be of the first 
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order, as the order parameters and excitations of the G phase in both spin and 
orbital sectors hardly change at such small doping levels. Specifically, AF spin cou- 
pling J c in the G phase may obtain a double-exchange ferromagnetic correction 
Jde < K/AS 2 = tx/2 due to the local vibrations of doped holes, which is about 
2 meV only at the critical doping x ~ 0.02. This correction is much smaller than 
J c ~ 6 me\S 

of undoped YVO3. Regarding the orbitals, they are perfectly aligned 
" ferromagnetically" along the c axis in the ground state of YVO3, satisfying the 
GdFeC>3 distortion; we find that this is not affected by 1-2% doping at all. Our 
predictions are then as follows: The staggered spin moment and spinwave disper- 
sions of the G phase remain barely unrenormalized upon doping, until the system 
suddenly changes its ground state. On the other hand, we expect that doping should 
lead to sizable changes in the properties of spin-C phase, because of the positive in- 
terplay between the "holons" formed and underlying orbital quantum physics. Our 
predictions are opposite to those of Ref.j^J this gives yet another opportunity to 
check experimentally (by means of neutron scattering) which approach — SE quan- 
tum picture or classical JT orbitals — is more appropriate for the spin-C phase of 
vanadates. 

A more intriguing question is, however, whether the doping induced C-phase of 
YVO3 is also dimerized as in the undoped YVO3 itself, or not. Theoretically, this is 
well possible at small doping, and, in fact, we expect that the orbital Peierls effect 
should cooperate with a conventional fermionic one, stemming from the hopping- 
modulation of doped carriers along the underlying orbital chains. 

Summarizing this section, we conclude that the ti g spin-orbital model with high 
spin values shows an intrinsic tendency towards dimerization. Considered together 
with lattice induced interactions, this leads to a several low-energy competing many- 
body states, and to the "fine" phase-selection by temperature, distortions and dop- 
ing. One of these phases is the spin-G phase with classical orbitals; it is well isolated 
in the Hilbert space and hence may show up only via a discontinuous transition. 

§5. Lifting the orbital degeneracy by spin-orbit coupling: Cobaltates 

It is well that in tig orbital systems a conventional spin-orbit cou- 

pling X(S ■ I) may sometimes play a crucial role in lifting the orbital degeneracy, 
particularly in case of late-transition metal ions with large A. When it dominates 
over the superexchange and weak orbital-lattice interactions, the spin and orbital de- 
grees of freedom are no longer separated, and it is more convenient to formulate the 
problem in terms of the total angular momentum. Quite often, the spin-orbit ground 
state has just the twofold Kramers degeneracy, and low-energy magnetic interactions 
can be mapped on the pseudospin one-half Hamiltonian. This greatly reduces the 
(initially) large spin-orbital Hilbert space. But, as there is "no free lunch", the pseu- 
dospin Hamiltonians obtain a nontrivial structure, because the bond directionality 
and frustrations of the orbital interactions are transfered to the pseudospin sector via 
the spin-orbital unification. Because of its composite nature, ordering of pseudospin 
necessarily implies both S- and Z-orderings, and the "ordered" orbital in this case 
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is a complex function. Via its orbital component, the pseudospin order pattern is 
rather sensitive to the lattice geometry. 

Experimental indications for such physics are as follows: (?) a separate JT- 
structural transition is suppressed but a large magnetostriction effects occur upon 
magnetic ordering; (ii) effective g— values in the ground state may deviate from the 
spin-only value and are anisotropic in general; (Hi) magnetic order may have a com- 
plicated nontrivial structure because of the non-Heisenberg symmetry of pseudospin 
interactions. As an example, the lowest level of Co 2+ ions (t\ g e^, S = 3/2, I = 1) in 
a canonical Mott insulators CoO and KC0F3 is well described by a pseudospin one- 
halfP Low-energy spin waves in this pseudospin sector, which are separated from 
less dispersive local transitions to the higher levels with different total momentum, 
have been observed in neutron scattering experiments in KCoFaPS 

In this section, we apply the pseudospin approach to study the magnetic corre- 
lations for the triangular lattice of C0O2 planes, in which the Co 4+ i° ns interact 
via the ~ 90° Co-O-Co bonds. We derive and analyze the exchange interactions 
in this system, and illustrate how a spin-orbit coupling leads to unusual magnetic 
correlations. This study is motivated by recent interest in layered cobalt oxides 
Na2;Co02 which have a complex phase diagram including superconductivity, charge 
and magnetic orderingsf^ 

The C0O2 layer consists of edge sharing CoOg octahedra slightly compressed 
along the trigonal (c || [111]) axis. Co ions form a 2D triangular lattice, sandwiched 
by oxygen layers. It is currently under debate,^ whether the undoped C0O2 layer 
(not yet studied experimentally) could be regarded as Mott insulator or not. Con- 
sidering the strongly correlated nature of the electronic states in doped cobaltates as 
is supported by numerous magnetic, transport and photoemission measurements, we 
assume that the undoped C0O2 plane is on the insulating side or near the borderline. 
Thus the notion of spin-charge energy scale separation and hence the superexchange 
picture is valid at least locally. 

5.1. Superexchange interaction for pseudospins 

A minimal model for cobaltates should include the orbital degeneracy of the 
Co 4+ — ion, ^3 where a hole in the d 5 (t2 g ) shell has the freedom to occupy one out 
of three orbitals o = d yz , b = d xz , c = d xy . The degeneracy is partially lifted by 
trigonal distortion, which stabilizes A\ g electronic state (a+b+c) / \/3 over ^-doublet 
{e ±i(t> a + e^b + c)/V3 (hereafter = 2tt/3): 

H A = A[n{E' g )-2n(A lg )]/3. (5-1) 

The value of A is not known; Reff^ estimates it A ~25 meV. Physically, A should 
be sample dependent thus being one of the control parameters. We will later see 
that magnetic correlations are indeed very sensitive to A. 

In terms of the effective angular momentum / = 1 of the t2g-shell, the functions 
Ai g and E' g correspond to the \l z = 0) and \l z = ±1) states, respectively. Therefore, 
a hole residing on the E' g orbital doublet will experience an unquenched spin-orbit 

interaction H\ = —\(S ■ I). The coupling constant A for a free Co 4+ ion is 650 
cm 1 ps 80meV :741 (this may change in a solid due to the covalency effects). 
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Fig. 4. (a) The i2 9 -orbital degeneracy of Co 4+ (d 5 )-ion is lifted by trigonal distortion and spin-orbit 
interaction. A hole with pseudospin one-half resides on the Kramers /-doublet. Its wavefunction 
contains both E' g and Ai g states mixed up by spin-orbit coupling, (b) Hopping geometry on 
the triangular lattice of Co-ions. a/3(7) on bonds should be read as t a ~ t, t~ n = —t', and 
a, 13, j £ {a,b,c} with a = d yz , b — d xz , c = d xy . The orbital nondiagonal t-hopping stems 
from the charge-transfer process via oxygen ions, while t' stands for the hopping between the 
same orbitals due to either their direct overlap or via two intermediate oxygen ions involving 
t pp . (After RefPS). 

The Hamiltonian H = Ha + H\ is diagonalized by the following transforma- 
tion!^ 

a a = {[cee-^'f-s + iseh + e^g 9 + sge^'h-s - ic e h a ]/V3 , (5-2) 

where cq = cos8,sg = sinO, a = (a,b,c) and ip a = (</>,— <f>,0), correspondingly. The 
angle 9 is determined from tan 26 = 2\/2A/(A + 2 A). As a result, one obtains three 
levels, fa, 9a, h 9 [see FigQJa)], each of them are Kramers doublets with pseudospin 
one-half a. The highest, /-level, which accommodates a hole in d? configuration, is 
separated from the (/-level by ej — e g = A + |(A/2 + A)(l/ cos 28—1). This splitting 
is ~ 3A/2 at A 3> A, and ~ A in the opposite limit. It is more convenient to use a 
hole representation, in which the level structure is reversed such that /-level is the 
lowest one. It is important to observe that the pseudospin f a states 

|f)/ = *c»| + l,i)-Sfl|0,T), 

|I)/ = *c 9 |-l,T)-sj|0 ) l> (5-3) 

are coherent mixture of different orbital and spin states. This will have important 
consequences for the symmetry of intersite interactions. 

We assume the hopping Hamiltonian suggested by the edge-shared structure: 

H? = t{a\j 3 a + 0la ja ) - t'llalja + h.c. , (5-4) 

where t = t\/A p d originates from d-p-d process via the charge-transfer gap A p d, and 
t' > is given either by direct d-d overlap or generated by indirect processes like 
d-p-p-d. On each bond, there are two orbitals active in the d-p-d process, while the 
third one is transfered via the i'-channel [Fig^b)]. For simplicity, we assume t' < t 
and neglect t' for a while. 
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There are three important superexchange processes: (a) Two holes meet each 
other at the Co-site; the excitation energy is U d - (b) Two holes meet each other at an 
intermediate oxygen-site; the excitation energy 2A pd + U p . (c) The oxygen electron 
is transferee! to an unoccupied e g shell and polarizes the spin of the t2 g level via 
the Hund's interaction, — 2Jjj(s e • st). This process is important because the e g — p 
hopping integral t a is larger than t n . The process (b), termed "a correlation effect" 
by Goodenoug h@ is expected to be stronger than contribution (a), as cobaltates 
belong to the charge-transfer insulators.^ 

These three virtual charge fluctuations give the following contributions: 



(a): 


A0i 


Sj + l/4)(n ia n jb + n ib n ja + ajftjatfy + feto^ta,-), 


(5-5) 


(b): 


B(Si 


Sj - l/A)(n ia n jb + n ib n ja ), 


(5-6) 


(c): 


- C(Si 


Sj){ n ic + n jc)> 


(5-7) 



where A = 4t 2 /U d , B = At 2 / (A pd +U p /2) and C = BR, with R ~ (2J H / A pd )(t a /U) 2 ; 
R ~ 1.5—2 might be a realistic estimation. [As usual, the orbitals involved depend on 
the bond direction, and the above equations refer to the 1-2 bond in FigEfb)]. The 
first, A— contribution can be presented in a SU(4) form Q3T5JI like in titanites and 
may have ferromagnetic as well as antiferromagnetic character in spin sector depend- 
ing on actual orbital correlations. While the second (third) contribution is definitely 
antiferromagnetic (ferromagnetic). As the constants A~i?~C~ 20 — 30 meV 
are smaller than spin-orbit splitting, we may now project the above superexchange 
Hamiltonian onto the lowest Kramers level /, obtaining the effective low-energy in- 
teractions between the pseudospins one-half Sf. Two limiting cases are presented 
and discussed below. 

Small trigonal field, A <C A. — The pseudospin Hamiltonian has the most sym- 
metric form when the quantization axes are along the Co-0 bonds. For the Co-Co 
pairs 1-2, 2-3 and 3-1 [Fig^b)], respectively, the interactions read as follows: 

H(l - 2) = J eff (-Sf Sf - SVS] + Sf Sf ), (5-8) 
H(2 - 3) = J eff (-S?Sy - StS z 3 + SfSf), 
H(3 - 1) = J ef f(-S?S? - SfS* + S*S?), 

where J e ff = 2(B + C)/9 ~ 10 — 15 meV. Here, we denoted the /-pseudospin simply 
by S. Interestingly, the ^4-term (|5-5|1 does not contribute to the /-level interactions 
in this limit. The projected exchange interaction is anisotropic and also strongly 
depends on the bond direction, reflecting the spin-orbital mixed character of the 
ground state wave function, see Ea. (|5-3|) . Alternation of the antiferromagnetic com- 
ponent from the bond to bond, superimposed on the frustrated nature of a triangular 
lattice, makes the Hamiltonian (|5-8[) rather nontrivial and interesting. Surprisingly, 
one can gauge away the minus signs in Eas. ([5-8|) in all the bonds simultaneously. 
To this end, we introduce four triangular sublattices, each having a doubled lattice 
parameter 2a. The first sublattice includes the origin, while the others are shifted 
by vectors 012, S23 and 031 [1,2,3 refer to the sites shown in FigHfb)]. Next, we 
introduce a fictitious spin on each sublattice (except for the first one), S', which are 



40 



G. Khaliullin 



obtained from the original /-pseudospins S by changing the sign of two appropri- 
ate components, depending on sublattice index [this is a similar trick, used in the 
context of the orbital Hamiltonian (|3-5|) for ferromagnetic YTiOa; see for details 
RefsJ^H23|]. After these transformations, we arrive at very simple result for the 
fictitious spins: H = Jeff 01 ■ S'j) in all the bonds. Thus, the known result 
for the AF-Heisenberg model on triangular lattice can be used. Therefore, we take 
120°— ordering pattern for the fictitious spins and map it back to the original spin 
space. The resulting magnetic order has a large unit cell shown in FigEI [Four sub- 
lattices have been introduced to map the model on a fictitious spin space; yet each 
sublattice contains three different orientations of S']. Ferro- and antiferromagnetic 
correlations are mixed up in this structure, and the first ones are more pronounced 
as expected from Ea. (|5-8|) . Magnetic order is highly noncollinear and also noncom- 
planar, and a condensation of the spin vorticity in the ground state is apparent. 
The corresponding Bragg peaks are obtained at positions K/2 = (2-7r/3, 0), that is, 
half-way from the ferromagnetic T— point to the AF-Heisenberg K— point. Because 
of the "hidden" symmetry (which becomes explicit and takes the form of a global 
SU(2) for the fictitious spins S"), the moments can be rotated at no energy cost. 
[Note that the "rotation rules" for the real moments are not that simple as for S': 
they are obtained from a global SU(2) rotation of S' via the mapping S <^=^ S']. 
Thus, the excitation spectrum is gapless (at A = 0), and can in fact be obtained by 
folding of the spinwave dispersions of the AF-Heisenberg model. We find nodes at 
the Bragg points K/2 (and also at r, K, M points, but magnetic excitations have a 
vanishing intensity at this points) . Spin- wave velocity is v = (3\/3/2) J e ff. A doped- 
hole motion in such spin background should be highly nontrivial — an interesting 
problem for a future study. 

Large trigonal field, A 3> A. — In this natural quantization axes are the 

in-plane (ab) and out-of-plane (c) directions. The low energy magnetic Hamiltonian 
obtained reads as follows: 

H(A » A) = J c SfSI + J ab (SfS* + Sj) , (5-9) 

with J c = [A - 2(3R - l)B]/9 and J ab = (A - B)/9. As A ~ B and R > 1, 
we have a large negative J c and small J ab (ferromagnetic Ising-like case). Thus, 
the large compression of C0O2 plane stabilizes a uniaxial ferromagnetic state with 
moments aligned along the c direction, and the excitations have a gap. When the 
AF coupling between the different planes is included, expected magnetic structure 
is then of A- type (ferro-Co02 planes coupled antiferromagnetically) . 

The above results show that the nature of magnetic correlations is highly sen- 
sitive to the trigonal distortion. This is because the variation of the ratio X/A 
strongly modifies the wave- function of the lowest Kramers doublet, which deter- 
mines the anisotropy of intersite exchange interactions. In this context, it is tempt- 
ing to mention recent NMR-reports on the observation of magnetic correlations in 
superconducting cobaltates at wave vectors "in between" the ferromagnetic and AF 
120-degree modesp*' while correlations change towards the ferromagnetic point in 
samples with a larger trigonal splitting. 801 These observations can naturally be un- 
derstood within the theory presented here, assuming that the Na-doped compounds 
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Fig. 5. Expected magnetic ordering in undoped C0O2 plane in the limit of a small trigonal field. 
Shown is the magnetic unit cell which contains twelve different lattice sites. Circles with ± sign 
indicate the out-of-plane component of the magnetic moment canted away from the plane by 
an angle 8 (with tan 8 = ±\/2). Associated Bragg spots in a momentum space are located at 
K/2 — (2tt/3, 0) and equivalent points. (No Bragg intensity at K points). Note that correlations 
in a majority of the bonds are more "ferromagnetic" rather than " antiferromagnetic" . 



still "remember" a local magnetic interactions that we obtained here for the insulat- 
ing limit. 

Pseudofermion pairing. — Another interesting point is that unusual superex- 
change interactions in C0O2 layer may have important consequences for the pairing 
symmetry in doped compounds, as argued recently in Ref. 751 Because of the non- 
Heisenberg form of the effective J— Hamiltonian, usual singlet/triplet classification 
is not very appropriate. To visualize a symmetry of two paired spins — in the 
spirit of the RVB picture — we represent the exchange Hamiltonian (|5-8j) in terms 
of /-fermionic spinons corresponding to the lowest Kramers doublet. Choosing now 
quantization along z || c, we find 

Ha = -J eff A^A%\ 4f = (t ijl0 + e^V + eT^VO/Vs . (5-10) 

Here, t ijfi = i(/;f /jj+Zij/jf )/\/2 , Uj,l = /jf/jf and % _i = ffj^ correspond to dif- 
ferent projections M = 0, ±1 of the total pseudospin Stot = 1 of the pair. The phase 
<t>^ in (T^TUl) depends on the (ij)-bond direction: = (0, 4>, ~4>) for (12, 23, 13)- 
bonds [see Fig|I[b)], respectively. As is evident from Eg. 1)5-10(1 . the pairing field 
is spin symmetric but nondegenerate, because it is made of a particular linear com- 
bination of M = 0, ±1 components, and the total spin of the pair is in fact quenched 
in the ground state. The absence of degeneracy is due to the fact that the pairing 
J— interaction has no rotational SU(2) symmetry. In a sense of spin non-degeneracy, 
the pair is in a singlet state, although its wavefunction is composed of M = 0,±1 
states and thus spin symmetric. [For the fictitious spins introduced above, this would 
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however read antisymmetric]. When such unusual pairs do condense, a momentum 
space wave- function must be of odd symmetry (of p, f, ..-wave character; a precise 
form of the gap functions is given in RefP). It is interesting to note, that the mag- 
netic susceptibility (and hence the NMR Knight shift) is partially suppressed in all 
three directions in this paired state. [We recall that S represents a total moment of 
the lowest Kramers level; but no Van Vleck contribution is considered here]. The rel- 
ative weights of different M-components, which control the Knight-shift anisotropy, 
depend on trigonal distortion via the ratio A/X, see for details RefP 

5.2. Spin/orbital polarons in a doped NaCo02 

Finally, we discuss an interesting consequence of the orbital degeneracy in sodium- 
rich compounds Nai_ x Co02 at small x. They are strongly correlated metals and also 
show magnetic order of ^4-type, which seems surprising in view of the fact that only 
a small number x of magnetic Co 4+ ions are present (in fact, magnetism disappears 
for large x). The parent compound NaCoC>2 is usually regarded as a band insulator, 
as Co 3+ ions have a spinless configuration pL. However, one should note that this 
state results from a delicate balance between the e g -t2 g crystal field splitting and a 
strong intraatomic Hund's interaction, and a control parameter lODq — 2Jh may 
change sign even under relatively weak perturbations. This would stabilize either 
^2g e g or ^2g e g magnetic configurations, bringing thereby Mott physics "hidden" in the 
ground state configuration. Doping of NaCo02 (by removing Na) is a very efficient 
way of doing it: a doped charge on Co + sites breaks locally the cubic symmetry 
and hence splits the e g and t2 g levels on neighboring Co 3+ sites. This reduces the 
gap between the lowest level of the split e g doublet and the upper level of the t2 g 
triplet, favoring a t\ g e g (S=l) configuration.^ As a result, a doped charge Co 4+ 
is dressed by a hexagon of spin-and-orbitally polarized Co 3+ ions, forming a local 
object which can be termed as a spin/orbital polaron. The idea of orbital polarons 
has been proposed and considered in detail in RefP*' in the context of weakly doped 
LaMnOs. The e g level splitting on sites next to a doped hole has been estimated to 
be as large as ~ 0.6 eV (in a perovskite structure with 180°-bonds), leading to a large 
binding energies and explaining the insulating nature of weakly doped ferromagnetic 
manganites. 

In fact, doping induced spin-state transmutations have been observed in per- 
ovskite LaCoC^pl in which a small amount of Sr impurities triggers magnetism. 
Because of their nontrivial magnetic response to the doping, we may classify NaCo02 
and LaCoOs as Mott insulators with incipient magnetism. Indeed, a ground state 
with a filled t2 g shell looks formally similar to that of band insulator but is qual- 
itatively different from the latter: NaCo02 and LaCo03 have low-lying magnetic 
states. In LaCoOs, they are just ~ 10 meV above the ground state. Thus, the spin 
and charge energy scales are completely different. In fact, Co 3+ ions in LaCo03 fully 
retain their atomic-like multiplet structure, as it is well documented by eseP and 
inelastic neutron scattering^ measurements. Spin states of such a Mott insulator 
can easily be activated by doping, temperature etc. 

The next important point is that the internal spin structure of the spin/orbital 
polarons in NaCo02 is very different from that in perovskites LaCo03 and LaMnC>3. 
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In the latter cases, polarons have a large spin due to internal motion of the bare hole 
within the polaron (local double-exchange process).^ We argue now that due to 90°- 
geometry of Co-O-Co bonds in NaCo02, the exchange interactions within a polaron 
are strongly antiferromagnetic, and thus a polaron has a total spin one-half only. 
(i) Consider first a superexchange between Co 3+ ions that surround a central Co 4+ 
and form a hexagon. They are coupled to each other via the 90°-superexchange. An 
antiferromagnetic interaction J(Si-Sj) between two neighboring Co 3+ spins (S=l) is 
mediated by virtual hoppings of electron between t2 g and e g orbitals, t = t a t n /A p d, 
and we find J ~ tP/E, where E is a relevant charge excitation energy (order of X1& 
or A p d + Up/2). Note that there will be also a weaker contribution from t2 g — t2 g 
hoppings (antiferromagnetic again) and some ferromagnetic corrections from the 
Hund's interaction (oc Jjj/U), which are neglected. Considering t ~ 0.2 eV and 
E ~ 3 eV, we estimate J could be as large as 10-15 meV. Below this energy scale, 
spins S=l of Co 3+ — hexagon form a singlet and are thus "hidden", but they do 
contribute to the high-temperature magnetic susceptibility in a form C/(T — 9) with 
a negative 9 ~ —J. (ii) Next, one may think that the double exchange process 
between the t2 g S=l/2 of the central Co 4+ and S=l of the Co 3+ would stabilize 
a large-spin polaron, but this is not the case. Surprisingly, this interaction is also 
of antiferromagnetic character, because there is no direct e g — e g transfer (which 
gives strong ferromagnetic polarons in perovskites with 180°-bonds like LaCo03 
and LaMnOs). Instead, there are two nearly equal double-exchange contributions 
which favor different spin alignment: (A) an e g electron of Co 3+ goes to (a singly 
occupied) t2 g level of Co 4+ . This process is possible when the Co 3+ -Co 4+ pair is 
antiferromagnetic. Another process is that (B) a t2 g electron of Co 3+ goes to an 
empty e g level of Co 4+ , which favors ferromagnetic bond via the Hund's coupling as 
for the usual double exchange. The hopping amplitude in both cases is just the same, 
t, but the kinetic energy gain is slightly larger for the antiferromagnetic configuration, 
because an extra energy ~ WDq (somewhat blocking a charge transfer) is involved 
in the second process. Therefore, the total double exchange within the pair Co 3+ - 
Co 4+ is weakly antiferromagnetic (a small fraction of t). We thus expect that each 
polaron brings about a free spin one-half in total, which contributes to the low-energy 
spin ordering. However, a polaron has an internal excitations to a multitude of its 
larger-spin states which, in principle, could be observed by neutron scattering as a 
nondispersive and broad magnetic modes (in addition to a propagating excitations 
in a spin one-half polaron sector). 

Our overall picture for Nai_ x Co02 at small x is that of a heavy spin/orbital 
polaron liquid. When two polarons overlap, their spins one-half are coupled fer- 
romagnetically via the mutual polarization of their spin clouds (Co 3+ — hexagons). 
Therefore, the polaron liquid may develop (itinerant) ferromagnetism within the 
C0O2 planes at low temperatures. The polarons may form clusters or charge-ordered 
patterns (near commensurate fillings) as discussed in Ref. 8 '-' However, the polaron 
picture breaks down at large doping x, and no magnetism is expected therefore in 
sodium-poor compounds. In other words, a heavy-fermion behavior and magnetism 
of weakly doped NaCo02 originate from the spin-state transmutation of Co 3+ ions 
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near doped holes, resulting in a narrow polaron bands. An apparent paradox of a 
large negative seen in susceptibility, and a small positive one inferred from spinwave 
dispersions 86 ' is a natural consequence of our picture: The former one reflects strong 
AF-couplings (active at large T) within the polaron, while the latter stems from a 
relatively weak ferromagnetic couplings between the spin-one-half polarons at low en- 
ergy limit, contributing to magnons. The in-plane coupling subtracted from magnons 
is found to be small (only of the order of the coupling between the planes) 
This is because it stems from a residual effective interactions between a dilute sys- 
tem of spin-one-half polarons. Our explanation implies that a quasi two-dimensional 
layered structure of cobaltates should show up in magnetic excitations of the magnet- 
ically dense compounds close to the undoped parent system C0O2. In that limit, we 
believe that \ J a b/J c \ S> 1, just like in a closely related, but regular magnet NaNiC>2 
with the same crystal and (A-type) magnetic structure, where \J a b/Jc\ ~ 15P^ An 
important remark concerning the magnon gaps: since all the orbital levels are well 
split within the polaron (lifting the degeneracy by orbital-charge coupling 821 ), the 
magnetic moment of the polaron is mostly of the spin origin. Thus, a low-energy 
magnetic excitations of a polaron liquid should be of the Heisenberg-type, as in fact 
observed in sodium-rich samplesP^'^ 

In our view, marked changes in the properties of cobaltates around Nao.7Co02 
are associated with a collapse of the polaron picture. In the magnetic sector, this 
implies a breakdown of dynamical separation on the strong-AF (weak-F) magnetic 
bonds within (in-between) the polarons, and a uniform distribution of the exchange 
interactions sets in. 

Summarizing results for the magnetic interactions in the two limiting cases con- 
sidered in this section: — a "pure" C0O2 and weakly doped NaCo02 — we conclude, 
that the nature of magnetic correlations in cobaltates is very rich and strongly sensi- 
tive to the composition. Description in terms of a simple AF-Heisenberg models with 
a uniform interactions in all the bonds is not sufficient. Depending on the strength 
of the trigonal distortion and doping level, we find the spin structures ranging from a 
highly nontrivial one as shown in Fig|SJ to a conventional A-type structures. The A- 
type correlations, which would be favored in a sodium-poor sample by large trigonal 
field, have an uniaxial anisotropy. On the other hand, the A-type state in a sodium- 
rich compounds is isotropic. Besides a conventional magnons, this state may reveal 
an interesting high-energy response stemming from its "internal" polaron structure. 

§6. Summary 

We considered several mechanisms lifting the orbital degeneracy, which are based 
on: (i) electron-lattice coupling, (it) the spin-orbital superexchange, and (Hi) a rel- 
ativists spin-orbit coupling. We discussed mostly limiting cases to understand spe- 
cific features of each mechanism. Reflecting the different nature of underlying forces, 
these mechanisms usually compete, and this may lead to nontrivial phase diagrams 
as discussed for example for YVO3. This underlines the important role of the or- 
bital degrees of freedom as a sensitive control parameter of the electronic phases in 
transition metal oxides. 
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We demonstrated the power of the "orbital-angle" ansatz in manganites, where 
the validity of classical orbital picture is indeed obvious because of large lattice 
distortions. In titanites and vanadates with much less distorted structures, we find 
that this simple ansatz is not sufficient and the orbitals have far more freedom in 
their "angle" dynamics. Here, the lattice distortions serve as an important control 
parameter and are thus essential, but the quantum exchange process between the 
orbitals becomes a "center of gravity". 

Comparing the e g and t<i g exchange models, we found that the former case is more 
classical: the e g orbitals are less frustrated and the order-from-disorder mechanism 
is very effective in lifting the frustration by opening the orbital gap. In fact, the 
low-energy fixed point, formed below the orbital gap in e g spin-orbital models, can 
qualitatively be well described in terms of classical "orbital- angle" picture. A strong 
JT nature of the e g quadrupole makes in reality this description just an ideal for all 
the relevant energy /temperature scales. 

A classical "orbital-angle" approach to the t2 g spin-orbital model fails in a fatal 
way. Here, the orbital frustration can only be resolved by strong quantum disorder 
effects, and dynamical coupling between the spins and orbitals is at the heart of the 
problem. This novel aspect of the "orbital physics", emerged from recent experimen- 
tal and theoretical studies of pseudocubic titanites, provide a key in understanding 
unusual properties of these materials, and make the field as interesting as the physics 
of quantum spin systems. We believe that the ideas developed in this paper should 
be relevant also in other ti g orbital systems, where the quantum orbital physics has 
a chance to survive against orbital-lattice coupling. 

Unusual magnetic orderings in layered cobaltates, stabilized by a relativistic 
spin-orbit coupling are predicted in this paper, illustrating the importance of this 
coupling in compounds based on late-3d and Ad ions. Finally, we considered how the 
orbital degeneracy can be lifted locally around doped holes, resulting in a formation 
of spin/orbital polarons in weakly doped NaCo02- The idea of a dilute polaron liquid 
provides here a coherent understanding of otherwise puzzling magnetic properties. 

In a broader context, a particular behavour of orbitals in the reference insulating 
compounds should have an important consequences for the orbital-related features of 
metal-insulator transitions under pressure/temperature/composition, e.g. whether 
they are " orbital-selective'^^ or not. In case of LaTiOs, where the three-band 
physics is well present already in the insulating state, it seems natural not to expect 
the orbital selection. On the other hand, the orbital selection by superexchange 
and/or lattice interactions is a common case in other insulators {e.g., xy orbital 
selection in vanadates); here, we believe that the "dimensionality reduction" phe- 
nomenon — which is so apparent in the "insulating" models we discussed — should 
show up already in a metallic side, partially lifting the orbital degeneracy and hence 
supporting the orbital-selective transition picture. 
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